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P , ABSTRACT 

We present a comprehensive analysis of weak gravitational lensing by large-scale structure in the Hubble Space Telescope Cosmic 
Evolution Survey (COSMOS), in which we combine space-based galaxy shape measurements with ground-based photometric red- 
shifts to study the redshift dependence of the lensing signal and constrain cosmological parameters. After applying our weak lensing- 
optimized data reduction, principal component interpolation for the spatially and temporally varying ACS point-spread function, 
' and improved modelling of charge-transfer inefficiency, we measure a lensing signal which is consistent with pure gravitational 

modes and no significant shape systematics. We carefully estimate the statistical uncertainty from simulated COSMOS-like fields 
, obtained from ray-tracing through the Millennium Simulation, including the full non-Gaussian sampling variance. We test our lens- 

ing pipeline on simulated space-based data, recalibrate non-linear power spectrum corrections using the ray-tracing analysis, employ 
photometric redshift information to reduce potential contamination by intrinsic galaxy alignments, and marginalize over system- 
atic uncertainties. We find that the weak lensing signal scales with redshift as expected from General Relativity for a concordance 
ACDM cosmology, including the full cross-correlations between different redshift bins. Assuming a flat ACDM cosmology, we mea- 
sure erg (Q. m /03) 03> = 0.75 ± 0.08 from lensing, in perfect agreement with WMAP-5, yielding joint constraints Q m = 0.266+qq23> 
cr 8 = 0.802+jJ JJ^ (all 68.3% conf.). Dropping the assumption of flatness and using priors from the HST Key Project and Big-Bang nu- 
cleosynthesis only, we find a negative deceleration parameter qo at 94.3% confidence from the tomographic lensing analysis, providing 
independent evidence for the accelerated expansion of the Universe. For a flat wCDM cosmology and prior w e [-2,0], we obtain 
w < -0.41 (90% conf). Our dark energy constraints are still relatively weak solely due to the limited area of COSMOS. However, 
' they provide an important demonstration for the usefulness of tomographic weak lensing measurements from space. 

H ; 
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1. Introduction 

During the last decade strong evidence for an acceler- 
ated expansion of the Universe has been found with sev- 
eral independent cosmological probes including type la su- 
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pernoyae ( Riess et al 



l2007t iKowalski et a" 



crowave background Jde Bernardis et al 



1998; Perlmutt eretaT]|1999l: iRiess et all 
120091), cosmic mi- 
2000 ; ISpergel et al l 
^slAllen etal.1 12008; 



20081: iHicken et a 



2003; Komatsu et al. 2009), galaxy clusters 
Mantz et alJ l2008[ 120091: IVikhlinin et alJ l2009h. barvon acous 
tic oscillations (Eisenstein et al 



2005; iPercival et al.1 



2009 ), integrated Sachs- Wolfe effect ( Giannant onio et al] 



2007 



2008 



Granett et all 120081: iHo et alJ 120081) . and strong gravitational 



lensing JSuvu et all 120091) . Within the standard cosmological 
framework this can be described with the ubiquitous presence 
of a new constituent named dark energy, which counteracts the 
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attractive force of gravity on the largest scales and contributes 
~ 70% to the total energy budget today. There are various at- 
tempts to explain dark energy, ranging from Einstein's cosmo- 
logical constant, via a dynamic fluid named quintessence, to a 
possible breakdown of Gen eral Relativity (e.g. lHuterer & Linderl 
2007i: lAlbrecht et al.ll20"09l) . all of which lead to profound impli- 
cations for fundamental physics. In order to make substantial 
progress and to be able to distinguish between the different sce- 
narios, several large dedicated surveys are currently being de- 
signed. 

One technique holding particularly high promise to con- 
strain dark energy jAlbrecht et alj l2"006t iPeacock eta 
Albr echt et al.l l2009h is weak gravitational lensing, which uti- 
lizes the subtle image distortions imposed onto the observed 
shapes of distant galaxies while their light bundles pass 
through the gravitational pote ntial of foreground structures (e.g. 
Bartelmann & Schneider 2001). The strength of the lensing ef- 
fect depends on the total foreground mass distribution, indepen- 
dent of the relative contributions of luminous and dark matter. 
Hence, it provides a unique tool to study the statistic al proper- 
ties of large-scale structure directly (for reviews see Sch neider! 
l2006tlHoekstra & Jairtl2008tlMunshi etafll2008h. 

Since its first detectio ns by Bac on et alJ (120001): iKaiser et al] 
(120001) ; IVan Waerbeke et al.l (120001) and IWittman et alf (l2000). 
substantial progress has been made with the measurement of 
this cosmological weak lensing effect, which is also called cos- 
mic shear. Larger sur veys have significan t ly reduced statisti- 
cal uncertainties (e.g. Hoekstra et al. 2002 1 iBrown et al .|l2003: 
Jarvis et alJ2"oMlMassev et al.l2005tlyan Waerbeke et al.l2005[ 
Hoekstra et al. 200fx ISemboloni et al.l l2006; Het terscheidt et al] 
2007; Fu et al. 2008), while tests on simu lated data have led to a 
better understanding of PSF systematic s dHevmans et al.ll2006at 
iMassev et al]|2007at iBridle et al.ll2010l and references therein). 
Finally, being a geometric effect, gravitational lensing depends 
on the source redshift distribution, where most earlier measure- 
ments had to rely on external redshift calibrations from the small 
Hubble Deep Fie lds. Here, the impact o f sampling variance was 
demonstrated by Benjam in et al.l who recalibrated ear- 

lier measurements using photometric redshifts from the much 
larger CFHTLS-Deep, significantly improving derived cosmo- 
logical constraints. 

Dark energy affects the distance-redshift relation and sup- 
presses the time-dependent growth of structures. Being sen- 
sitive to both effects, weak lensing is a powerful probe of 
dark energy properties, also prov iding important tests for the- 
ories of modified gravity (e.g. Benabed & Bernard eaul 1200 it 



Benabed & van Wae rbeke 2004; Schi md et alJl2007t iDore et al. 
20071; Uain & Zhangj 120081; ISchmidtl 120081) . Yet, in order to 
significantly constrain these redshift-dependent effects, the 
shear signal must be measured as a function of source 
redshift, an analysis often calle d weak lensing t omogra- 



20021: iJain & Taylor! I2003T iHeavensI 12003b Hu & Jainl 



phy o r 3D weak lensing (e.g. iHul 1 19991 12001 IMu terer 



2004 



Bernstein & Jainll2004l; ISimon et alj|2004t iTakada & Jainll2004 
Heavens et al .1 12006: Ta ylor et al.l 120071) . Redshift information 



is additionally required to eliminate potential contamination 
of the lensing signal f r om intrinsic galaxy alignments (e.g. 
King & Schneider) 120021: iHirata & Seliakl l2004t IHevmans etall 
2006bj iJoachimi & Schneider! 120081) . In general, weak lensing 
studies have to rely on photometric reds hifts (e.g. lBenite3 [2000; 
lllbert et al.ll2006HHildebrandt et al.ll2008l) given that most of the 
studied galaxies are too faint for spectroscopic measurements. 

So far, tomographic c osmological weak l e nsing techniques 
were applied to real data bv lBacon et al.ld2005l) : ISemboloni et all 



(120061) : iKitching et all (120071) : IMassev et al] (l2007d) . Dark en- 
ergy constraints from previous weak lensing surveys were lim- 
ited by the lack of the required individual photometric redshifts 
dJarvis et al.ll2006t iHoekstra et al.ll2006l: ISemboloni et alJl2 006: 
Kilbi nger et alJ l2009al) or small survey area ( Kitching et al.l 
2007). The currently best data- set for 3D weak l ensing is 



given by the COSMOS Survey dScoville et all l2007al) . which 
is the largest continuous area ever imaged with the Hubble 
Space Telescope (HST), comprising 1.64 deg 2 of deep imag- 
ing with the Advanced Camera for Surveys (ACS). Compared 
to ground-based measurements, the HST point-spread func- 
tion (PSF) yields substantially increased number densities of 
sufficiently resolved galaxies and better control for systemat- 
ics due to smaller PSF corrections. Although HST has been 
used for earlier co s mological weak lensing analyses (e.g . 



[ Refregier et al. 20021 iRhodes et all l2004t iMiralles et al.l 12005b 
IHevmans et al. 120051; ISchrabback et al J 120071) . these studies 
lacked the area and deep photometric redshifts which are avail- 
able for COSMOS (lllbert et alJl2009l) . This combination of su- 
perb space-based imaging and ground-based photometric red- 
shifts make s COSMOS the perfe ct test case for 3D weak lens- 
ing studies. IMassev et"a"D d2007d) conducted an earlier 3D weak 
lensing analysis of COSMOS, in which they correlated the shear 
signal between three redshift bins and constrained the matter 
density Q m and power spectrum normalization cr 8 . In this paper 
we present a new analysis of the data, with several differences 
compared to the earlier study: we employ a new, exposure-based 
model for the spatially and temporally varying ACS PSF, which 
has been derived from dense stellar fields using a principal com- 
ponent analysis (PCA). Our new parametric correction for the 
impact of charge transfer inefficiency (CTI) on stellar images 
eliminates earlier PSF modelling uncertainties caused by confu- 
sion of CTI- and PSF-induced stellar ellipti city. Using the latest 
photometric redshift catalogue of the field (lllbert et al. 2009), we 
split our galaxy sample into five individual redshift bins and ad- 
ditionally estimate the redshift distribution for very faint galax- 
ies forming a sixth bin without individual photometric redshifts, 
doubling the number of galaxies used in our cosmological anal- 
ysis. We study the redshift scaling of the shear signal between 
these six bins in detail, employ an accurate covariance matrix 
obtained from ray-tracing through the Millennium Simulation, 
which we also use to recalibrate non-linear power spectrum cor- 
rections, and marginalize over parameter uncertainties. In addi- 
tion to Q m and <x 8 , we also constrain the dark energy equation of 
state parameter w for a flat wCDM cosmology, and the vacuum 
energy density £2a for a general (non-flat) ACDM cosmology, 
yielding constraints for the deceleration parameter qo. 

This paper is organized as follows. We summarize the most 
important information on the data and photometric redshift cat- 
alogue in Sect.|2l while further details on the ACS data reduc- 
tion are given in App.lAl Section[3]summarizes the weak lensing 
measurements including our new correction schemes for PSF 
and CTI, for which we provide details in App.|B] We conduct 
various tests for shear-related systematics in Sect.|4] We then 
present the weak lensing tomography analysis in Sect. [5] and 
cosmological parameter estimation in Sect.|6] We discuss our 
findings and conclude in Sect. [7] 

Throughout this paper all magnitudes ar e given in the AB 
system , where z'gw denotes the SExtractor dBertin & Arnoutsl 
1996) MAGJVUTO magnitude measured from the ACS data 
(Sect.12. 11>. while i + is the MAG_AUT0 magnitude determined by 
lllbert et al J (120091) from the Subaru data (S ect.l2!2~TT) In severa l 
tests we employ a reference WMAP-5-like (Dunkley et al. 2009) 
flat ACDM cosmology characterized by Q m = 0.25, <x 8 = 0.8, 
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h = 0.7 2, Qb = 0-044, n s — 0.9 6, where we use the transfer func- 
tion by Eisenst ein & Hul(|l998l) and non-li near power spectrum 
corrections according to Smit h et al. (200j). 

2. Data 

2.1. HST/ACS data 

The COSMOS Survey dScoville et al.ll2007al) is the largest con- 
tiguous field observed with the Hubble Space Telescope, span- 
ning a total area of ~ 77' x 77' (1.64deg 2 ). It comprises 579 
ACS tiles, each observed in F814W for 2028s using four 
dithered exposures. The survey is centred at a = 10 h 00 m 28.6 s , 
6 = +02°12'21'/0 (J2000.0), and data were taken between 
October 2003 and November 2005. 

We have reduced the ACS/WFC data starting from the 
flat-fielded images. We apply updated bad pixel masks, sub- 
tract the sky background, and compute optimal weights as de- 
tailed in App.fA] For the image registration, distortion correc- 
tion, cosmic ray rejectio n, and stacking we use MultiDrizzle^ 
dKoekemoer et al.l2002|). applying the latest time-dependent dis- 
tortion solution from I Anderson! (l2007l) . We iteratively align ex- 
posures within each tile by cross-correlating the positions of 
compact sources and applying residual shifts and rotations. 

In tests with dense stellar fields we found that the default 
cosmic ray rejection parameters of MultiDrizzle can lead to 
false flagging of central stellar pixels as cosmic rays, especially 
if telescope breathing introduces significant PSF variations (see 
Sect.[3jl between combined exposures. Hence, stars will be par- 
tially rejected in exposures with deviating PSF properties. On 
the contrary, galaxies will not be flagged due to their shallower 
light profiles, leading to different effective stacked PSFs for 
stars and galaxies. To avoid any influence on the lensing anal- 
ysis, we create separate stacks for the shape measurement of 
galaxies and stars, where we use close to default cosmic ray 
rejection parameters for the for mer (driz_cr_snr="4 . 3.0" , 
driz_cr_scale="1.2 0.7", see iKoekemoer et all 120021 I2007I) . 
but less aggressive masking for the latter (driz_cr_snr="5.0 
3.0", driz_cr_scale="3.0 0.7"). As a result, the false masking 
of stars is substantially reduced. On the downside some actual 
cosmic rays lead to imperfectly corrected artifacts in the "stel- 
lar" stacks. This is not problematic given the very low fraction 
of affected stars, for which the artifacts only introduce additional 
noise in the shape measurement. 

For the final image stacking we employ the LANCZ0S3 in- 
terpolation kernel and a pixel scale of 0705, which minimizes 
noise correlations and aliasing without unnecessarily broaden- 
ing the PSF ( for a detailed comparison to other kernels see 
iJee et alj [2007). Based on our input noise models (see App.fA} 
we compute a correctly scaled RMS image for the stack. We 
match the s t acked image WCS to the ground-based catalogue by 
lllbert et alj d2009l) . 

We employ our RMS noise model for o bject detec- 
tion with SExtractor dBertin & Arnoutsl I1996I) . where we 
require a minimum of 8 adjacent pixels being at least 
1.4cr above the background, employ deblending parameters 
DEBLENDJJTHRESH = 32, DEBLENDJIINCONT = 0.01, and mea- 
sure MAGJUJTO magnitudes 1^4, which we correct fo r a mean 
galactic extinction offset of 0.035 (iSchlegel et al.ll998l) . Objects 
near the field boundaries or containing noisy pixels, for which 
fewer than two good input exposures contribute, are auto- 
matically excluded. We also create magnitude-scaled polyg- 
onal masks for saturated stars and their diffraction spikes. 

1 MultiDrizzle version 3.1.0 



Furthermore, we reject scattered light and large, potentially in- 
correctly deblended galaxies by running SExtractor with a low 
0.5cr detection threshold for 3960 adjacent pixels, where we fur- 
ther expand each object mask by six pixels. The combined masks 
for the stacks were visually inspected and adapted if necessary. 

Our fully filtered mosaic shear catalogue contains a to- 
tal of 446 934 galaxies with z'gi4 < 26.7, corresponding to 
76 galaxies/arcmin 2 , where we exclude double detections in 
overlapping tiles and reject the fainter component in the case 
of close galaxy pairs with separations < 075. For details on the 
weak lensing galaxy selection criteria see App. IB.6l 

In addition to the stacked images, our fully time-dependent 
PSF analysis (see Sect. [3] App. IB.5t makes use of individual ex- 
posures, for which we use the cosmic ray-cleansed COR images 
before resampling, provided by MultiDrizzle during the run 
with less aggressive cosmic ray masking. These are only used 
for the analysis of high signal-to-noise stars, which can be iden- 
tified automatically in the half-light radius versus signal-to-noise 
spacfl Here we employ simplified field masks only excluding 
the outer regions of a tile with poor cosmic ray masking. 

2.2. Photometric redshifts 

2.2.1. Individual photometric redshifts for i + < 25 galaxies 

We use th e public COSMO S-30 photometric redshift cata- 
logue from lllbert et alJ d2009h . which covers the full ACS mo- 
saic and is magnitude limited to z' + < 25 (Subaru SExtractor 
MAG_AUT0 magnitude). It is based on the 30 band photomet- 
ric catalogue, which includes imaging in 20 optical bands, as 
well as near- infrared and deep IRAC data (Capak et al. 2009, in 
preparation), lllbert et alj d2009l) computed photomet ric redshift 
using the he Phare code (S. Arnouts & O. Ilbert; also lllbert et al.l 
2006), reaching an excellent accuracy of cr A ,/( 1+Z ) = 0.012 for 
i + < 24 and z < 1.25. The near-infrared (NIR) and infrared cov- 
erage extends the capability for reliable photo-z estimation to 
higher redshifts, where the B a imer break moves out of the op- 
tical bands. Extended to z ~ 2. lllbert etal] (12009) find an accu- 
racy of cr A ,/ (I+; . s ) = 0.06 at z' + ~ 24. The comparison to spectro- 
scopi c redshifts from the zCOSMOS-deep sample (Lill y et all 
120071) with z'jL dian = 23.8 indicates a 20% catastrophic out- 
lier rate (defined as |z p hot -Zs P ed/(l +z sp ec) > 0.15) for galax- 
ies at 1.5 < Zspec < 3. In particular, for 7% of the high-redshift 
fapec > 1 -5) galaxies a low-redshift photo-z (z p hot < 0.5) was as- 
signed. This degeneracy is expected for faint (z' + > 24) high- 
redshift galaxies, for which the Balmer break cannot be iden- 
tified if they are undetected in the NIR data (limiting depth 
J ~ 23.7, K ~ 23.7at5cr). Due to the employed magnitude prior 
the contamination is expected to be mostly uni-directional from 
high to low redshifts. 

We tested this by comparing the COSMOS-30 ca t alogue 
to photometric redshifts estimated by iHildebrandt et al.1 (2009) 
in the overlapping CFHTLS-D2 field using onl y optical u*gri z 
bands and the BPZ photometric redshift code (Benitez 2000). 
Here we indeed find that 56% of the matched z' + < 25 galaxies 
with COSMOS-30 photo-zs in the range 2 < zc30 < 4 are identi- 
fied at zd2 < 0.6 in the D2 catalogue, if only a weak cut to reject 
galaxies with double-peaked D2 photo-z PDFs (ODDS > 0.7) is 
applied^- 



2 Ar h = 0.25 pixel wide kernel; S /N > 20, defined as in Erben et al. 
d2001h : peak flux < 25 000e~ 

3 A more stringent cut ODDS > 0.95 reduces this fraction to 14%. 
Yet, it also reduces the absolute number of galaxies by a factor 4.7. Note 
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If not accounted for, such a contamination of a low-photo-z 
sample with high-redshift galaxies would be particularly severe 
for weak lensing tomography, given the strong dependence of the 
lensing signal on redshift. In Sect.[5]we will therefore split galax- 
ies with assigned z p hot < 0.6 into sub-samples with expected low 
(z' + < 24) and high (z' + > 24) contamination, where we only in- 
clude the former in the cosmological analysis. Matching our 
shear catalogue to the fully masked COSMOS-30 photo-z cat- 
alogue yields a total of 194 976 unique matches. 

2.2.2. Estimating the redshift distribution for z' + > 25 galaxies 

In order to include galaxies without individual photo-zs in our 
analysis, we need to estimate their redshift distribution. Fig.Q] 
shows the mean photometric COSMOS-30 redshift for galaxies 
in our shear catalogue as a function of z'gi4. In the whole magni- 
tude range 23 < igu < 25 the data are very well described by the 
relation 



<z> = (0.276 + 0.003)(/ 8 i4 - 23) + 0.762 + 0.003 



(1) 



For comparison we also plot points from the H ubble Deep Field- 
North (HDF-N, iFernandez-Soto et alj[l999l) and Hubble Ultra 
Deep Field (HUDF lCoe et alJ l2006fl for the extended magni- 
tude range 23 < i%u < 27, where both catalogues are redshift 
complete. The HDF-N data agree very well with the COSMOS 
fit over the whole extended range, on average to 2%. In con- 
trast, the mean photometric redshifts in the HUDF are on av- 
erage higher than ([]} by 16% for 23 < i%u < 25 and 10% for 
25 < z' 8 i4 < 27. The difference between the HDF-N and HUDF 
can be regarded as a rough estimate for the impact of sampling 
variance in such small fields. The fact that the HUDF galaxies 
systematically deviate from (Q]) not only for z'gi4 > 25 but also 
for /gi4 < 25 where COSMOS-30 photo-zs are available, indi- 
cates that it is most likely affected by sampling variance con- 
taining a relative galaxy over-density at higher redshift. Given 
the excellent fit for the COSMOS galaxies and very good agree- 
ment for the HDF-N data we are thus confident to use (Q]) for a 
limited extrapolation to < 26.7 for our shear galaxies. This 
is also motivated by the fact that igi4 < 25 and r'gi4 > 25 galax- 
ies are not completely independent, but partially probe the same 
large-scale structure at different luminosities. 

Due to the non-linear dependence of the shear signal on red- 
shift it is not only necessary to estimate the correct mean redshift 
of the galaxies, but also their actual redshift distribution. In weak 
lensing studies the redshift distribution is often parametrized as 
p(z) oc (z/zp) 0, exp [ -(z/zq)^] (e.g. iBrainerd et all 19961) . which 
ISchrabback et alJ y007) extended by fitting a,/3 in combina- 
tion with a linear dependence of the median redshift on mag- 
nitude, leading to a magnitude-dependent zo- Yet, it was noted 
that this fit was not fully capable to reproduce the shape of the 
redshift distribution of the fitted galaxies. Given the higher ac- 
curacy needed for the analysis of the larger COSMOS data, we 
use a modified parametrization 



p{z\ku) °c — 

\Z0 



exp 



+ cu exp 



(2) 



that, in contrast, 26% (22%) of the matched galaxies with a D2 photo-z 
2 < zd2 < 4 are placed at z C 30 < 0.6 for ODDS > 0.7 (ODDS > 0.95). 
These could be explained by Lyman-break galaxies, which are better 
constrained by the deeper u* observations in the CFHTLS-D2. In any 
case we expect negligible influence on our results given our treatment 
for faint zc30 ^ 0.6 galaxies. 

4 For the H UDF we interpolat e j 814 from the i lls and z 850 magnitudes 
provided in the lCoe et al catalogue. 
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Fig. 1. Relation between the mean photometric redshift and £§14 
magnitude for COSMOS, HUDF, and HDF-N, where the error- 
bars indicate the error of the mean assuming Gaussian scat- 
ter and neglecting sampling variance. The best fit ([T} to the 
COSMOS data from z'gi4 < 25 is shown as the bold line, whereas 
the thin lines indicate the conservative 10% uncertainty con- 
sidered for the extrapolation in the cosmological analysis. The 
HDF-N data agree with the relation very well, whereas the 
mean redshifts are higher in the HUDF both for z'gi4 < 25 and 
z'gi4 > 25, demonstrating the influence of sampling variance in 
such small fields. 



where zo = ZoO'814), and u = max[0, (z'su - 23)]. Using a max- 
imum likelihood fi0 we determine best-fitting parameters 
(a,/3,c,d,y) = (0.678,5.606,0.581, 1.851, 1.464) from the indi- 
vidual magnitudes, photo-zs, and (symmetric) 68% photo-z er- 
rors of all galaxies with 23 < z'gn < 25. From Eqs. (Q]) and (|2]) 
we then numerically compute the non-linear relation between zo 
and z'gi4, for which we provide the fitting formulae 

ZO = 0.446(z 8 i4- 23)+ 1.235 for 22 < ; 814 < 23 (3) 

7=7 

zo = ~ 23 )/ 4 ]' for 23 < '814 < 27 (4) 

7=0 

with (a ,...,a 7 ) = (1.237,1.691,-12.167,43.591,-76.076, 
72.567,-35.959,7.289). The total redshift distribution of 
the survey is then simply given by the mean distribution 
<f>(z) = 'LtiP(z\ i su,k)/N. 

We chose the functional form of © because its first addend 
allows for a good description of the peak of the redshift distribu- 
tion, while the second addend fits the magnitude-dependent tail 
at higher redshifts; see Fig.|2]for a comparison of the data and 
model in four magnitude bins. In Fig. [3] we compare the actual 
redshift distribution for the combined HDF-N and HUDF data 
to the one we predict from their magnitude distribution and the 
fit to the COSMOS data, finding very good agreement also for 
25 < z'gi 4 < 27. The only major deviation is given by a galaxy 
over-density in the HUDF photo-z data near z ~ 3.2, which is 



5 We employ the CERN Program Library 
( http : //wwwasdoc . web . cern . ch/wwwasdoc/minuit/ 1. 



MINUIT 



Tim Schrabback et al.: Evidence for the accelerated expansion of the Universe from 3D weak lensing with COSMOS 



5 




24 < i + < 25 are truly at z = 2, increasing the estimated mean 
redshift by 8%. Compared to the fit uncertainty in (Q]i (~ 1%) 
this constitutes the main source of error for our redshift extrap- 
olation. In the cosmological parameter estimation (Sect.[6]l, we 
constrain this uncertainty and marginalize over it using a nui- 
sance parameter, which rescales the redshift distribution within 
a conservatively chosen +10% interval. Note that the +4% dif- 
ference between the measured and predicted mean redshift of the 
combined HDF-N and HUDF data in Fig. [3] actually suggests a 
smaller uncertainty. 




Fig. 2. Redshift histogram for galaxies in our shear catalogue 
with COSMOS-30 photo-zs (dotted), split into four magnitude 
bins. The solid curves show the fit according to (0 and (|2j, 
which is capable to describe both the peak and high redshift tail. 




Fig. 3. Combined redshift histogram for the HDF-N and HUDF 
photo-zs, split into two magnitude bins. The solid curves show 
the prediction according to ((TJ, (0 and the galaxy magnitude dis- 
tribution. The good agreement for 25 < igu < 27 galaxies con- 
firms the applicability of the model in this magnitude regime. 



also partially responsible for the higher mean redshift in Fig.Q] 
and which may be attributed to large-scale structure. 

Our fitting scheme assumes that the COSMOS-30 photo-zs 
provide unbiased estimates for the true galaxy redshifts. 
However, in Sect. 12.21 we suspected that z' + > 24 galaxies with 
assigned z < 0.6 might contain a significant contamination with 
high-redshift galaxies. To assess the impact of this uncertainty, 
we derive the fits for ([TJ and (|2]l using only galaxies with 
23 < i + < 24, reducing the estimated mean redshift of shear 
galaxies without COSMOS-30 photo-z by 4%. As an alterna- 
tive test, we assume that 20% of the z < 0.6 galaxies with 



3. Weak lensing shape measurements 

To measure an accurate lensing signal, we have to carefully cor- 
rect for instrumental signatures. Even with the high-resolution 
space-based data at hand, we have to accurately account for both 
PSF blurring and ellipticity, which introduce spurious shape dis- 
tortions. To do so, one requires both a good model for the PSF, 
and a method which accurately employs it to measure unbi- 
ased estimates for the (reduced) gravitational shear g from noisy 

galaxy images. 

Fo r the latter, we use the K S B+ formalism dKaiser et al.l 
1 19951: iLuppino & Kaiserl Il997t iHoekstra et alJ 1 19981) . see 
Erbe n et al.1 d2001l) : ISchrabback et al.l (120071) and App.|BJ] for 
details on our implementation. As found with simulations of 
ground-based weak lensing data, KS B+ can significantly un- 
deres t imate gravitational she ar (Erben et al. 20 0U iBacon et al.l 
2001; iHevmans et al.l l2006at iMassev et al.l l2007al) . where the 
calibration bias m and possible PSF anisotropy residuals c, de- 
fined via 



^obs 



(5) 



depend on the details of the implementation. IMassev et al.l 
(2007a, STEP2) detected a shear measurement degradation for 
faint objects for our pipeline, which is not surprising given the 
fact that the KSB+ form alism does not account for noise. While 
ISchra bback et al. (2007) simply corrected for the resulting mean 
calibration bias, the 3D weak lensing analysis performed here 
requires unbiased shape measurements not only on average, but 
also as f unction of redshift, and hence galaxy magnitude and size 
(see e.g. lKitohingetal] l2008. 20091 ISemboloni et alJl2009h . We 
therefore empirically account for this degradation with a power- 
law fit to the signal-to-noise dependence of the calibration bias 



m = -0.078 



S/N 
2 



-0.38 



(6) 



where S/N is c omputed with the galaxy size-dependent KSB 
weight function (lErben et al.ll200ll), and co rrected for noise cor- 
relations as done in lHartlap et al.l (12009). As S/N relates to 
the significance of the galaxy shape measurement, it provides 
a more direct correction for noise-related bias than fits as a 
function of magnitude or size. We have determined this correc- 
tion usin g the STEP2 simula tions of ground-based weak lens- 
ing data (Massev et al. 2007a). In order to test if it performs re- 
liably for the ACS data, we have analysed a set of simulated 
ACS-like data (see App. IB.2l i. In summary, we find that the re- 
maining calibration bias is m = +0.008 + 0.002 on average, and 
\m\ < 0.02 over the entire magnitude range used, which is neg- 
ligible compared to the statistical uncertainty for COSMOS. 
Likewise, PSF anisotropy residuals, which are characterized in 
<(3j by c, are found to be negligible in the simulation (dispersion 
cr c = 0.0006), assuming accurate PSF interpolation. 
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Weak lensing analyses usually create PSF models from the 
observed images of stars, which have to be interpolated for the 
position of each galaxy. Typically, a high galactic latitude ACS 
field contains only ~ 10-20 stars with sufficient S /N, which are 
too few for the spatial polynomial interpolation commonly used 
in ground-based weak lensing studies. In addition, a stable PSF 
model cannot be used, given that substantial temporal PSF vari- 
ations have been detected, mostly caused by focus changes re- 
sulting from orbital temperature variations (telescope breathing), 
mid-term seasonal effects, and lon g-term shrinkage of the opti- 
cal telescope assembly (OTA) (e.g. Krist 2003; Lallo et al. 2006; 
Anderson & King! 12005 ISchrabback et al] l2007t iRhodes etal] 
2007). To circumvent this problem, we have implemented a 
PSF correction scheme based on principal compon ent analy- 
sis (PCA), as first suggested by iJarvis & Jainl (|2004). We have 
analysed 700 i%u exposures of dense stellar fields, interpolated 
the PSF variation in each exposure with polynomials, and per- 
formed a PCA analysis of the polynomial coefficient variation. 
We find that ~ 97% of the total PSF ellipticity variation in 
random pointings can be described with a single parameter re- 
lated tojhechajigejntelescope focus, confirming earlier results 
(e.g. lRhodes et al. 2007). However, we find that additional vari- 
ations are still significant. In particular, we detect a dependence 
on the relative angle between the pointing and the orbital tele- 
scope movement! suggesting that heating in the sunlight does 
not only change the telescope focus, but also creates slight addi- 
tional aberrations dependent on the relative sun angle. These de- 
viations may be coherent between COSMOS tiles observed un- 
der similar orbital conditions. To account for this effect, we split 
the COSMOS data into 24 epochs of observations taken closely 
in time, and determine a low-order, focus-dependent residual 
model from all stars within one epoch. We provide further de- 
tails on our PSF correction scheme in App. |B.5l 

As an additional observational challenge, the COSMOS data 
suffer from defects in the ACS CCDs, which are caused by the 
continuous cosmic ray bombardment in space. These defects act 
as charge traps reducing the charge-transfer-efficiency (CTE), an 
effect referred to as charge-transfer-inefficiency (CTI). When the 
image of an object is transferred across such a defect during par- 
allel read-out, a fraction of its charge is trapped and statistically 
released, effectively creatin g charge-trails following objects in 
the read-out y-direction (e .g. lRhodes et al.l2007t lChiaberee et al. 
2009; Massey et al.ll20Toh . For weak lensing measurements the 
dominant effect of CTI is the introduction of a spurious ellip- 
ticity component in the read-out direction. In contrast to PSF 
effects, CTI affects objects non-linearly due to the limited depth 
of charge traps. H ence, the two effects must be corrected sepa- 
rately. As done by Rho des et al.l |2007), we employ an empiri- 
cal correction for galaxy shapes, but also take the dependence 
on sky background into account. Making use of the CTI flux- 
dependence, we additionally determine and apply a paramet- 
ric CTI model for stars, which is important as PSF and CTI- 
induced ellipticity get mixed otherwise. We present details on 
our CTI correction sche mes for stars i n App.lB.4l and for galaxies 
in App. |B.6l Note that iMassev et al.l (|2010) recently presented 
a method to correct for CTI directly on the image level. We 
find that the methods employed here are sufficient for our sci- 
ence analysis, as also confirmed by the tests presented in Sect.|4] 
However, for weak lensing data with much stronger CTE degra- 
dation, such as ACS data taken after Servicing Mission 4, their 
pixel-based correction should be superior. 



4. 2D shear-shear correlations and tests for 
systematics 

To measure the cosmological signal and conduct tests for sys- 
tematics we compute the second-order shear-shear correlations 



U0) = 



(7) 



from galaxy pairs separated by § — \§i — Here, A (J - = 1 if 
the galaxy separation § falls within the considered angular 
bin around 9, and Ay = otherwise. In (01 we approximate 
our reduced shear estimates g = y/(l — k) ^ y with the shear 
y as commonly done in cosmological weak lensing (typically 
\k\ ~ 1% - 3%; correction employed in Sect. 16.41 ). decompose it 
into the tangential component y t and the 45 degree rotated cross- 
component y x relatively to the separation vector, and employ 
uniform weights. 

As an important consistency check in weak gravitational 
lensing, the signal can be decomposed into a curl-free compo- 
nent (E-mode) and a curl component (B-mode). Given that lens- 
ing creates only E-modes, the detection of a significant B-mode 
indic ates the presence of unc orrected residual systematics in the 
data. Crittenden et al.l ((2002) show that £ ± can be decomposed 
into E- and B-modes as 



M0) = 

with 



Jo ~» 



(8) 



(9) 



We plot this decomposition for our COSMOS catalogue in the 
left panel of Fig. [4] Given that the integration in © extends to 
infinity, we employ ACDM predictions for 6 > 40', leading to 
a slight model-dependence, which is indicated by the dashed 
curves corresponding to cr% = (0.7, 0.9), whereas the points have 
been computed for cr 8 = 0.8. Within this section, error-bars and 
covariances are estimated from 300 bootstrap resamples of our 
galaxy shear catalogue, which accounts for both shot noise and 
shape noise. As seen in Fig.|4] we detect no significant B-mode 
£;b- However, note that different angular scales are highly corre- 
lated for^£/B, which mixes power on a broad range of scales and 
potentially smears out the signatures of systematics. 

An E/B-mode decomposition, for which the correlation be- 
tween different sc ales is weaker, is provided by the dispersion of 
the aperture mass (ISchneiderl Q 996 ) 



i r 26 



d#4? 
IF 



f + (#)r + |^|±f_(#)r_|~ 



(10) 



6 Technically speaking, we show a dependence on the velocity aber- 
ration plate scale factor in Fig. IB.6l 



with T ± given in Sch neider et al.l (|2002), wh ere we employ 
t he ap erture mass weight function proposed bv lSchneider et al] 
d!998l) . The computation of ( fTUb requires integration from 
zero, which is not practical for real data. We therefore 
truncate £ ± for 6 < 0f05, where t he introduced bias i s 
small compar e d to ou r statistical errors dKilbinger et al.ll2006h . 
Mass ev et alJ d2007cl) measure a significant B-mode compo- 
nent {M 2 ) at scales 1' < 6 < 3', whereas this signal is neg- 
ligible in the present analysis. We quantify the on average 
slightly positive {M\) by fitting a mean offset taking the boot- 
strap covariance into account (correlation between neighbour- 
ing points =2 0.5), yielding an average {M 2 ) — (1 + 4) x 10~ 7 if 
all points are considered, and {M]){6 < 6') = (1.0 ± 1.4) X 10~ 6 




Fig. 4. Decomposition of the shear field into E- and B-modes using the shear correlation function £e/b (left), aperture mass dispersion 
(M^Ljj) (middle), and ring statistics CR!R.)e/b (right). Error-bars have been computed from 300 bootstrap resamples of the shear 
catalogue, accounting for shape and shot noise, but not for sampling variance. The solid curves indicate model predictions for 
erg = (0.7, 0.8). In all cases the B-mode is consistent with zero, confirming the success of our correction for instrumental effects. 
For %e/b the E/B-mode decomposition is model-dependent, where we have assumed <x 8 = 0.8 for the points, while the dashed 
curves have been computed for <x 8 = (0.7, 0.9). The dotted curves indicate the signal if the residual ellipticity correction discussed 
in App. IB.6l is not applied, yielding nearly unchanged results. Note that the correlation between points is strongest for £e/b and 
weakest for CRR) E /b- 



or (M])(8 < 2') = (4.0 + 4.7) x 10~ 6 if only small scales are in- 
cluded, consistent with no B-modes. 

The c leanest E/B-mode decomposi t ion is given by the ring 
statis tics dSchneider & Kilbingeri 120071: lEifler et al.ll2009bt see 
also [Fu & Kilbingeri 1201 Ol) . which can be computed from the 
correlation function using a finite interval with non-zero lower 
integration limit 



1 f*d& 

2 X* T [ ^ 



(#)Z + (#,77)±£_(0)Z_(#,77)] , (11) 



with functions Z ± given in ISchneider & Kilbingeri d2007l) . We 
compute CRf^ F./B using a scale -dependent integration limit 77 as 
outlined in Ei fler et alj d2009bl) . As can be seen from the right 
panel of Fig.|4j also CRf^s is consistent with no B-mode signal. 

The non-detection of significant B-modes in our shear cat- 
alogue is an important confirmation for our correction schemes 
for instrumental effects and suggests that the measured signal is 
truly of cosmological origin. 

As a final test for shear-related systematics we compute the 
correlation between corrected galaxy shear estimates y and un- 
corrected stellar ellipticities e* 



(r,/x< /x >Krr/x< /x )l 



/X^f/X 



(12) 



which we normalize using t he stellar auto-correlation as sug- 
gested by Bacon et al.l ([2003). As detailed in App. IB.6l we em- 
ploy a somewhat ad hoc residual correction for a very weak 
remaining instrumental signal. We find that £ sys is indeed only 
consistent with zero if this correction is applied (Fig. [5]), yet 
even without correction, £ sys is negligible compared to the ex- 
pected cosmological signal. The negligible impact can also be 
seen from the two-point statistics in Fig. [4] where the points are 
computed including residual correction, while the dotted lines 
indicate the measurement without it. We suspect that this resid- 
ual instrumental signature could either be caused by the limited 
capability of KSB+ to fully correct for a complex space-based 
PSF, or a residual PSF modelling uncertainty due to the low 
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Fig. 5. Cross-correlation between galaxy shear estimates and un- 
corrected stellar ellipticities as defined in (fT2l . The signal is con- 
sistent with zero if the residual ellipticity correction discussed in 
App. IB.6l is applied (circles). Even without this correction (trian- 
gles) it is at a level negligible compared to the expected cosmo- 
logical signal (dotted curves), except for the largest scales, where 
the error-budget is anyway dominated by sampling variance. 



number of stars per ACS field. In any case we have verified that 
this residual correction has a negligible impact on the cosmolog- 
ical parameter estimation in Sect.|6] changing our constraints on 
<x 8 at the 2% level, well within the statistical uncertainty. 
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Table 1. Definition of redshift bins, number of contributing 
galaxies, and mean redshifts. 



Bin 


Zmin 


'max 


N 


<z> 


1 


0.0 


0.6 


i + < 24 : 22294* 


0.37 








i + > 24 : 29 817 




2 


0.6 


1.0 


58 194 


0.80 


3 


1.0 


1.3 


36 382 


1.16 


4 


1.3 


2.0 


25 928 


1.60 


5 


2.0 


4.0 


21718 


2.61 


6 


0.0 


5.0 


251958 


1.54 ±0.15 



*: Here we also exclude 259 galaxies with ( + < 24, which have a 
significant secondary peak in their redshift probability distribution at 
Zphot.2 > 0.6. 



5. Weak lensing tomography 

In this section we present our analysis of the redshift dependence 
of the lensing signal in COSMOS. We start with the definition 
of redshift bins in Sect. 15. II summarize the theoretical frame- 
work in Sect. 15. 21 describe our angular binning and treatment of 
intrinsic galaxy alignments in Sect. 15. 31 elaborate on the covari- 
ance estimation in Sect. 15.41 present the measured redshift scal- 
ing in Sect. 15.51 and discuss indications for a contamination of 
faint Zphot < 0.6 galaxies with high redshift galaxies in Sect.| 



5.1. Redshift binning 

We split the galaxies with individual COSMOS-30 photo-zs into 
five redshift bins, as summarized in Table [TJ and illustrated in 
Fig. [6] We chose the intermediate limits z = (0.6, 1.0, 1.3) such 
that the Balmer/4000A break is approximately located at the 
centre of one of the broadband r + i + z + filters. This minimizes 
the impact of possible artifical clustering in photo-z space and 
hence scatter between redshift bins for galaxies too faint to be 
detected in the Subaru medium bands. Given our chosen limits, 
most catastrophic redshift errors are faint bin 5 galaxies identi- 
fied as bin 1 (~Sect. l2.27TT >. Thus, we do not include z < 0.6 galax- 
ies with i + > 24 in our analysis due to their potential contami- 
nation with high redshift galaxies, but study their lensing signal 
separately in Sect. 15. 61 We use all galaxies without individual 
photo-z estimates with 22 < i&u < 26. "fl as a broad bin 6, for 
which we estimated the redshift distribution in Sect. 12.2.21 

5.2. Theoretical description 

Extending the formalism from Sect. [4] we split the galaxy sam- 
ple into redshift bins and cross-correlate shear estimates between 
bins k and / 



X,,Ay 



(13) 



where the summation extends over all galaxies i in bin k, and 
all galaxies j in bin /. These are estimates for the shear cross- 
correlation functions which are filtered versions of the con- 
vergence cross-power spectra 



i r 

2* Jo 



(14) 



7 Including galaxies with i + < 25 which are located in masked re- 
gions for the ground-based photo-z catalogue, but not for the space- 
based lensing catalogue. 




Fig. 6. Redshift distributions for our tomography analysis. The 
solid-line histogram shows the individual COSMOS-30 red- 
shifts used for bins 1 to 5, while the difference between the 
dashed and solid histograms indicates the 24 < z' + < 25 galax- 
ies with Zphot < 0.6, which are excluded in our analysis due to 
potential contamination with high-redshift galaxies. The long- 
dashed curve corresponds to the estimated redshift distribution 
for 

!8i4 < 26.7 shear galaxies without individual COSMOS-30 
photo-z, which we use as bin 6. 

where J„ denotes the n th -order Bessel function of the first kind 
and I is the modulus of the two-dimensional wave vector. These 
can be calculated from line-of-sight integrals over the three- 
dimensional (non-linear) power spectrum Pg (see Sect. l6.2b as 



()"m 



4c 4 



X 



d * i77\ — Pg 

a z (x) 
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X 



(15) 



with the Hubble parameter Hq, matter density D. m , scale factor a, 
comoving radial distance^-, comoving distance to the horizon;^, 
and comoving angular diameter distance /k(x)- The geometric 
lens-efficiency factors 



gk(x) 



Jk(x'-x) 



4*r' Pkix')- , . 
x Mx) 



are 
the 



weighted according to the redshift distributions 
two considered redshift bins (see e.g. 



(16) 



of 



iBartelmann & Schne ider 200lT: ISimon et al.| [2004). 
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5.3. Angular binning and treatment of intrinsic galaxy 
alignments 

Our six redshift bins define a total of 21 combinations of redshift 
bin pairs (including auto-correlations). For each redshift bin pair 
(k, I), we compute the shear cross-correlations £f and in six 
logarithmic angular bins between 0'2 and 30'. We include all 
of these angular and redshift bin combinations in the analysis 
of the weak lensing redshift scaling presented in this section, to 
keep it as general as possible. Yet, for the cosmological parame- 
ter estimation in Sect. [6] we carefully select the included bins to 
minimize potential bias by intrinsic galaxy alignments and un- 
certainties in theoretical model predictions. 
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In order to minimize potential contamination by intrinsic 
alignments of physically associated galaxies, we exclude the 
auto-correlations of the relatively narrow redshift bins 1 to 5. 
These contain the highest fraction of galaxy pairs at similar red- 
shift, and hence carry the strongest potential contamination. 

An additional contamination may originate from alignments 
between intrinsic galaxy shapes and the ir surrounding density 
field causing the g ravitational shear (e.g. IHirata & Sel iak 2004; 
iHirata et al.ll2007l) . A complete removal of this effect requires 
more advanced analysis schemes (e.g. Ijoachimi & Schneider] 
2008), which we postpone to a future st udy. Yet, following the 
suggestion by Mandelba um et all (120061) . we exclude luminous 
red galaxies (LRGs) in the computation of the shear-shear cor- 
relations used for the parameter estimation. This reduces po- 
tential contamination, given t hat LRGs were found to carry 
the stronges t alignment signal (Mandelbaum et alj |2006i 120091; 
IHirata et al]|2007l) . We select these galaxies from the lllbert et al.l 
(2009) photo-z catalogue with cuts in the photometric type 
modgai < 8 ("ellipticals") and absolute magnitude My < -19, 
excluding a total of 5 75 1 galaxies^. We accordingly adapt the 
redshift distribution for the parameter estimation. 

In the cosmological parameter estimation, we additionally 
exclude the smallest angular bin (0 < 0.'5), for which the the- 
oretical model predictions have the largest uncertainty due to 
required non -linear correction s (Sect. l6.2i i and the influence of 
baryons (e.g. Rudd et al. 20Q8|). 

While we do not exclude LRGs and the smallest angular bin 
for the redshift scaling analysis presented in the current section, 
we have verified that their exclusion leads to only very small 
changes, which are well within the statistical errors and do not 
affect our conclusions. 

5.4. Covahance estimation 

In order to interpret our measurement and constrain cosmo- 
logical parameters, we need to re liably estimat e the da ta co- 
variance matrix and its inverse. iMassev et al.1 (l2007d) esti- 
mate a covariance for their analysis from the variation be- 
tween the four COSMOS quadrants. This approach yields too 
few independent real isations and may s ubstantially underesti- 
mate the true errors (Hartla p et al.ll2007h . We a lso do not em- 
ploy a covariance for Gaussian statistics (e.g. iJoachimi et al.l 
2008) due to the neglected influence of non-Gaussian sam- 
pling variance. This is particularl y important for the small-scale 
signal probed with C OSMOS dKilbinger & Schneider! 12001 
Sem boloni et al.ll2007l) . Instead, we estimate the covariance ma- 
trix from 288 realisations of COSMOS-like fields obtained from 
ray-tr acing through the Millennium Simulation ( Springel et al. 
2005), which combines a large simulated volume yielding many 
quasi-independent lines-of-sight with a relatively high spatial 
and mass resolution. The latter is needed to fully utilize the 
small-scale signal measureable in a deep space-based survey. 

The details o f the ray-tracing analysis are given in 
Hilb ertet al.l d2009). In brief, we use tilted lines-of-sight through 
the simulation to avoid repetition of structures along the 
backwards lightcone, providing us with 32 quasi-independent 
4 deg x4 deg fields, which we further subdivide into nine 
COSMOS-like subfields, yielding a total of 288 realisations. We 
randomly populate the fields with galaxies, employing the same 
galaxy number density, field masks, shape noise, and redshift 

8 In the cross-correlation between two redshift bins, it would be suf- 
ficient to exclude LRGs in the lower redshift bin only. However, for 
convenience we generally exclude them. 
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distribution as in the COSMOS data. We incorporate photomet- 
ric redshift errors for bins 1 to 5 by randomly misplacing galaxy 
redshifts assuming a (symmetric) Gaussian scatter according to 
the lcr errors in the photo-z catalogue. In contrast, the redshift 
calibration uncertainty for bin 6 is not a stochastic but a sys- 
tematic error, which we account for in the cosmological model 
fitting in Sect.|6] 

The value of cr% = 0.9 used for the Millennium Simulation 
is slightly high compared to current estimates. This will lead to 
an overestimation of the errors, hence our analysis can be con- 
sidered slightly conservative. We have to neglect the cosmology 
dependence of the covariance (Eifler et al. 2009a) in the parame- 
ter estimation, given that we have currently only one simulation 
with high resolution and large volume at hand. 

We need to invert the covariance matrix for the cosmological 
parameter estimation in Sect. [6] While the covariance estimate 
C» from the ray-tracing realizations is unbiased, a bias is intro- 
duced by correlated noise in the matrix inversion. To obtain an 
unbiased estimate for the inverse covariance C , we apply the 
correction 

C" 1 =cC.~' = n ~ P ~ 2 C*~ l foxp<n-2 (17) 
n — 1 

discussed in Hartla p et al.l (12007), where n — 288 is the number 
of independent realisations and p is the dimension of the data 
vector. As discussed in Sect. 15.31 we exclude the smallest angular 
bin and auto-correlations of redshift bins 1 to 5, yielding p = 160 
and a moderate correction factor c 0.4390. In contrast, for the 
full data vector including all bins and correlations (p = 252), a 
very substantial correction factor c ^ 0.1 185 would be required. 
Hence, our optimized data vector also leads to a more robust 
covariance inversion. 

In order to limit the required correction for the covariance 
inversion, we do not include more angular bins in our analy- 
sis. We have therefore optimized the bin limits using Gaussian 
covariances dJoachimi et al.l l2008) and a Fisher-matrix analysis 
aiming at maximal sensitivity to cosmological parameters. 

5.5. Redshift scaling of shear-shear cross-correlations 

We plot the shear-shear cross-correlations £* 6 between all red- 
shift bins and the broad bin 6 in Fig. [7] These cross-correlations 
carry the lowest shot noise and shape noise due to the large num- 
ber of galaxies in bin 6. The good agreement between the data 
and ACDM model already indicates that the weak lensing sig- 
nal roughly scales with redshift as expected. The errors corre- 
spond to the square root of the diagonal elements of the full ray- 
tracing covariance. Points are correlated not only within a red- 
shift bin pair, but also between different redshift combinations, 
as their lensing signal is partially caused by the same foreground 
structures. In addition, galaxies in bin 6 contribute to different 
cross-correlations. Note that our relatively broad angular bins 
lead to a significant variation of the theoretical models within a 
bin. When computing an average model prediction for a bin, we 
therefore weight according to the ^-dependent number of galaxy 
pairs within this bin. Likewise, we plot points at their effective 6, 
which has been weighted accordingly. 

Instead of plotting 21 separation-dependent, noisy cross- 
correlations, we condense the information into a single plot 
showing the redshift dependence of the signal. Here we assume 
that the predictions for our reference cosmology describe the rel- 
ative angular dependence of the signal sufficiently well, and fit 
the data points as 

^,flt ( ) = ^, re i^,mod (0)) (lg) 
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Fig. 7. Shear-shear cross-correlations £f between bins 1 to 6 
and bin 6, where points are plotted at their effective 8, weighted 
within one bin according to the ^-dependent number of con- 
tributing galaxy pairs. The curves indicate ACDM predictions 
for our reference cosmology with erg = 0.8. Corresponding 
points and curves have been equally offset along the x-axis for 
clarity. The error-bars correspond to the square root of the di- 
agonal elements of the full ray-tracing covariance. Note that the 
points are substantially correlated both between angular and red- 
shift bins, leading to the smaller scatter than naively expected 
from the error-bars. 



where £ k J' mod (8) is the model for the refer ence cosmology with 
erg - 0.8, and £ ± ' le is the fitted relative amplitude. In this fit, 
we take the full ray-tracing covariance between the angular 
scales into account. We plot the resulting 21 "collapsed" cross- 
correlations for both £ + and in Fig. [8] as a function of their 
model prediction at a reference angular scale of 0f8, where 
points are again correlated. For both cases the redshift scaling 
of the signal is fully consistent with ACDM expectations, show- 
ing a strong increase with redshift. This demonstrates that 3D 
weak lensing does indeed perform as expected. We note that 
for the signal is somewhat low for lower redshift combi- 
nations (smaller ^'< mod ) 5 whereas it is slightly increased com- 
pared to predictions at higher redshifts. This behaviour is not 
surprising as mo st massive structures in COSMOS are located 
at 0.7 < z < 0.9 dScoville et al.ll2007bl) . which create a lensing 
signal only for the higher redshift source bins. Slight differences 
between £ + and are also expected, given that they probe the 
power spectrum with different filter functions, see Eq. (1141 . 



5.6. Contamination of the excluded faint z < 0.6 sample with 
high-z galaxies 



As discussed in Sect. 12. 21 we expect a significant fraction of faint 
i + > 24 galaxies with assigned photometric redshift z p hot < 0.6 
to be truly located at high redshifts z true > 2. To test this hypoth- 
esis, we plot the collapsed shear cross-correlations for differ- 
ent samples of galaxies with assigned z p hot < 0.6 in Fig. [9] For 
the i + < 24 galaxies used in the cosmological analysis the sig- 
nal is well consistent with expectations, suggesting negligible 



contamination. For a 24 < z' + < 25 sample with single-peaked 
photo-z probability distribution a mild increase is detected. This 
is still consistent with expectations, suggesting at most low con- 
tamination. We also study a sample of galaxies each of which 
has a significant secondary peak in their photometric redshift 
probability distribution at z p hot,2 > 0.6, amounting to 36% of 
all 24 < z' + < 25 galaxies with z p h t < 0.6. This sample shows a 
strong boost in the lensing signal, suggesting strong contamina- 
tion with high-redshift galaxies. 

We can obtain a rough estimate for this contamination if we 
assume that the shear signal does actually scale as in our ref- 
erence ACDM cosmology. For simplicity we assume that the 
cross-contamination can be described as a uni-directional scatter 
from bin 5 to bin 1, and that the true redshifts of the misplaced 
galaxies follow the distribution within bin 5. The expected con- 
taminated signal is then given as a linear superposition of the 
cross-correlation predictions with bin 1 and bin 5 respectively, 
according to the relative number of contributing galaxy pairs 



j-ll,cont /i \2/-ll,mod , 2j-55,mod , r» /i 



15, mod 



(19) 



z-l/,cont / 1 \ >-l/,mod , >-/5,mod 

£+ = C 1 - r K+ + >f + 



for I > 1 , 



where r is the contamination fraction, i.e. the fraction of the bin 
1 galaxies with 24 < i + < 25 and a significant secondary peak 
in their photo-z PDF, which should have been placed into bin 
5. We fit the measured shear-shear cross-correlations with 
( fT9l as a function of r, where we fix the reference ACDM 
cosmology and employ a special ray-tracing covariance (gen- 
erated for r = 0.5), yielding an estimate for the contamination 
r = 0.7 + 0.2 (stat.) + 0.1 (sys.), where the systematic error indi- 
cates the response to a change in cr 8 by 0. 1 . This translates to a 
total contamination of (25 ± 7 + 4)% for the 24 < z' + < 25 galax- 
ies with z p hot < 0.6, which is consistent with our estimate for the 
redshift calibration uncertainty for bin 6 (Sect. l2.2^2l . Note that 
we also measure an increased signal in f_ for the sample with 
secondary photometric redshift peak, but do not include it in the 
fit ( fT9l due to the stronger deviations for £ w,rel in Fig. [8] An ad- 
equate inclusion would then require a more complex analysis 
scheme, with a comparison not to the model predictions, but to 
all measured cross-correlations. 

Our analysis provides an in teresting confirma tion for the 
photometric redshift analysis by Ilbe rt et al] d2009l) . which ap- 
parently succeeds in identifying sub-samples of (mostly) un- 
contaminated and potentially contaminated galaxies quite effi- 
ciently. 



6. Constraints on cosmological parameters 

6.1. Parameter estimation and considered cosmological 
models 

The statistical analysis of the shear tomography correlation 
functions, assembled as data vector d, is based on a standard 
Bayesian approach (e.g. MacKav 2003). Therein, prior knowl- 
edge of model parameters p is combined with the information 
on those parameters inferred from the new observation and ex- 
pressed as posterior probability distribution function (PDF) of 
P- 



P(p\d) = 



P(d\p)P(p) 
P(d) 



(20) 



Here, P(p) is the prior based on theoretical constraints and pre- 
vious observations, and P(d) denotes the evidence. The likeli- 
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Fig. 8. Shear-shear redshift scaling for £ + (left) and (right). Each point corresponds to one redshift bin combination, where 
we have combined different angular scales by fitting the signal amplitude ^* /lel relative to the model prediction £± ,mod (ff) for our 
reference ACDM cosmology with cr% = 0.8. The lower plots show the relative amplitude as a function of the model prediction 
£*'' mod (0!8) for a reference angular bin centred at — 0'8, whereas the amplitude has been scaled with ^' mod (0'.8) for the upper 
plots. Symbols of one kind correspond to cross-correlations of one bin with all higher-numbered bins. Within one symbol the partner 
redshift bins sort according to the mean lensing efficiency, from left to right as 1,2, 3, 6, 4, 5. Note that points are correlated as each 
redshift bin is used for six bin combinations, and given that foreground structures contribute to the signal of all bin combinations 
at higher redshift. The error-bars are computed from the full ray-tracing covariance, accounting for this influence of large-scale 
structure. 



hood function P(d\p) is the statistical model of the measurement 
noise, for which we choose a Gaussian model 

\nP(d\p) = -i [d- m(p)] t C- 1 [d - m(p)] + const, (21) 

where m(p) is the parameter-dependent model, and C~ l the in- 
verse covariance, which we estimated from the ray-tracing real- 
izations in Sect. 15. 41 

In our analysis we consider different cosmologi- 
cal models, which are characterized by the parameters 
p = (£2de> Qm> c 8s h, w, fz), with the dark energy density Que, 
matter density £2 m , power spectrum normalization o~%, Hubble 
parameter h, and (constant) dark energy equation of state 
parameter w. Here, /, denotes a nuisance parameter encap- 
sulating the uncertainty in the redshift calibration for bin 6 
as pe(z, fz) = p6(fzZ), which was discussed in Sect. 12.2.21 We 
consider 

- a flat ACDM cosmology with fixed w — - 1 , Q m e [0,1], and 

^DE — — 1 _ Qm, 

- a general (non-flat) ACDM cosmology with fixed w — -1 
and Q. DE = Q A e [0, 2], Q m e [0, 1.6], and 

- a flat wCDM cosmology with w e [-2,0], Q. m e [0, 1], and 

In all cases, we employ priors with flat PDFs for <r 8 e [0.2, 1.5] 
and f z 6 [0.9, 1.1]. In our default analysis scheme we also apply 
a Gaussian prior for h - 0.72 + 0.025, and assume a fixed baryon 
dens ity Qb _ 0.044 and spectral index n s = 0.96 as consistent 
with lDunklev et all d2009l) . where the small uncertainties on Ob 



and « s are negligible for our analysis. Note that we relax these 
priors for parts of the analysis in Sect. l6.3.2l and Sect. 16.41 



The practical challenge of the parameter estimation is to 
evaluate the posterior within a reasonable time, as the com- 
putation of one model vector for shear tomography correla- 
tions is time-intensive. For an efficient sampling of the pa- 
rameter space, we empl oy the Popul a tion M onte Carlo (PMC) 
method as described in Wra ith et al.1 d2009l). This algorithm i s 
an adaptive importance-sampling technique dCappe et alJl2007l) : 
instead of creating a sample under the posterior as done in tra- 
ditional Monte-Carlo M arkov chain (MCMC) techniques (e.g. 
IChris tensen et ^Lll200lb . points are sampled from a simple dis- 
tribution, the so-called proposal, in our case a mixture of eight 
Gaussians. Each point is then weighted by the ratio of the pro- 
posal to the posterior at that point. In a number of iterative steps, 
the proposal function is adapted to give better and better approx- 
imations to the posterior. We run the PMC algorithm for up to 
eight iterations, using 5 000 sample points in each iteration. To 
reduce the Monte-Carlo variance, we use larger samples with 
10000 to 20000 points for the final iteration. These are used 
to create density histograms, mean parameter values, and confi- 
dence regions. Depending on the experiment, the effective sam- 
ple size of the final importance sample was between 7 500 and 
17 700. We also cross-checked parts of the analysis with an inde- 
pendently developed code which is based on the traditional but 
less efficient MCMC approach, finding fully consistent results. 
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Fig. 9. Shear-shear redshift scaling for £J; as in Fig. [8] but now 
only cross-correlations with bin 1 (z < 0.6) are shown, hence 
the different axis scale. The signal from the z' + < 24 galaxies 
used in our cosmological analysis (crosses), is well consistent 
with the ACDM prediction (curve). Galaxies with 24 < i + < 25 
and a single-peaked photo-z probability distribution (circles) 
show a mildly increased but still consistent signal. In con- 
trast, 24 < i + < 25 galaxies with a significant secondary peak at 
Zphot,2 > 0.6 in their individual photo-z probability distribution, 
show a strong signal excess (squares), suggesting strong con- 
tamination with high-redshift galaxies. 



6.2. Non-linear power spectrum corrections 

To calculate model predictions for the correlation functions ac- 
cording to (fT4l i. ([TBT l. and (fTSI l. we need to evaluate the in- 
volved distance ratios and compute the non-linear power spec- 
trum Pg{k,z). Given a set of parameter values, the computa- 
tion of the distances and the linearly extrapolated power spec- 
trum is straightforward . We employ the transfer function by 
Eisen stein & Hul ([1998) for the latter, taking baryon damping but 
no oscillations into account ('shape fit'). 

For ACDM models we estimate t he ful l non-linear power 
spectru m according to ISmith et al] (120031) . [McDonald et al.l 
(2006) also provide non-linear power spectrum corrections 
for w + -1, but these were tested for a narrow range in 
erg = 0.897 + 0.097 only. We want to keep our analysis as gen- 
eral as possible, not having to assume such a strong prior on <x 8 . 
Following the icosmo code (Refregier et al. 2008) we i nstead 
interpolate the non-linear corrections from ISmith et aTl {2003) 
between the cases of a ACDM cosmology (w = -1) and an 
OCDM cosmology, acting as a dark energy with w = -1/3. This 
is achieved by replacing the parameter / = Qa /(1 - Qm) in the 
halo model fitting function (Smith et al] |2003h . This parameter 
is used to interpolate between spatially flat models with dark en- 
ergy (/ = 1) and an open Universe without dark energy (/ = 0). 
We substitute f by a new parameter f = -0.5(3w+ 1). Thus, we 
obtain /' = 1 for ACDM and /' = for wCDM with w = -1/3, 
mimicking an OCDM cosmology for which the original param- 
eter / vanished as well. 

To test this simplistic approximation, we compare the com- 
puted corrections for w = (-0.5,-1.5) to the fitting formulae 
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Fig. 10. Comparison of the fit formulae for the non-linear 
growth of structure in wCDM cosmologies. Shown is the three- 
dimensional matter power spectrum, normalized by the corre- 
sponding ACDM power spectrum, as a function of the wave 
vector k. In the upper panel we consider a wCDM cosmology 
with w = -0.5, in the lower panel o ne with w = -1.5. Soli d 
curves show the fit to the simulations bv lMcDonald et ail (12006). 
while the dashed l ines have been obtained by interpolating the 
Smith et al. (2003) fitting formulae between the cases of an 
OCDM and a ACDM cosmology as outlined in Sect. 16.21 Each 
fit formula has been computed at redshifts z — (black), z - 0.5 
(blue), and z = 1 (orange). While deviations are substantial at 
z = 0, the lensing analysis of the deep COSMOS data is mostly 
sensitive to structures at z > 0.4, where deviations are reason- 
ably small. Note that the remaining cosmological parameters 
have been set to their default WMAP5-like values, except for 
<r% = 0.9. 



from iMcDonald et all (120061) in Fig.[l0] Note that we use our 
fiducial cosmological parameters to obtain the se curves, except 
for cr s = 0.9, to match cr 8 = 0.897 ± 0.097 from lMcDonald et al.l 
(2006). For most of the scales probed by our measurement the 
two descriptions agree reasonably well. The modification of the 
halo fit follows the fits to the simulations more accurately on 
large scales and at higher redshif t, while it does not reproduce 
the tendency of the fits bv lMcDonald et all ((2006) to drop off for 
large wave vectors. The precision of the modification outlined 
above is sufficient for our aim to provide a proof of concept for 
weak lensing dark energy measurements. However, future mea- 
surements with larger data sets will require accurate fitting for- 
mulae for general w cosmologies. 

6.3. Cosmological constraints from COSMOS 
6.3.1. Flat ACDM cosmology 

We plot our constraints on Q m and <x 8 for a flat ACDM cos- 
mology and our default 3D lensing analysis scheme in Fig.fTTI 
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Table 2. Constraints on erg (Q m /0.3) ff , Q m , Qde, and w from the COSMOS data for different cosmological models and analysis 
schemes, using our default priors. We quote the marginali zed mean and 68.3 % confidence limits (16th and 84th percentiles) as- 
suming non-linear power spectrum corrections according to lSmith et alj d2003l) and the description given in Sect. 16.21 Our analysis 
of the Millennium Simulation (Sect. 16.4] ) suggests that the cr 8 -estimates should be reduced by a factor x0.95 due to biased model 
predictions for the non-linear power spectrum and reduced shear corrections. The power-law slopes a have typical fit uncertainties 
ofcr a -0.02. 



Cosmology Analysis a erg (n ra /0.3) ff Cl m CIde w 

Flat ACDM 3D 051 0.79 ± 0.09 0.32+^* 068^ ^ =1 

Flat ACDM 2D 0.62 0.68 ±0.11 030*™ a70 -044 _1 

General ACDM 3D 0.77 0.74 ±0.12 0.43™™ - 97 -qI _1 

Flat wCDM 3D 0.47 0.79 ± 0.09 0.30+° : ff O7oj^ -1.23^g 




£2m 

Fig. 11. Comparison of our constraints on £2 m and cr 8 for a flat 
ACDM cosmology using a 3D (blue solid contours) versus a 
2D weak lensing analysis (green dashed contours). The contours 
show the 68.3% and 95.4% credibility regions, where we have 
marginalized over the parameters which are not shown. The 2D 
analysis favours slightly lower cr 8 resulting from the lack of mas- 
sive structures in the field at low redshifts. Nonetheless, the con- 
straints are fully consistent as our ray-tracing covariance prop- 
erly accounts for sampling variance. 

(solid contours), showing the typical 'banana-shaped' degener- 
acy, from which we compuuf] 

cr 8 (n m /0.3) 51 = 0.79 + 0.09 (68.3% conf.). 

Here we marginalize over the uncertainties in h and the pa- 
rameter f, encapsulating the uncertainty in the redshift calibra- 
tion for bin 6, where we find that f z is nearly uncorrected with 
£l m , and only weakly correlated with cr 8 . The data allow us to 
weakly constrain f. = 1.03^ ^, with a maximum posterior point 
at /. = 1.05. This constraint is nearly unchanged for the other 
cosmological models considered below. 

For comparison we also conduct a classic 2D lensing anal- 
ysis (dashed contours in Fig.fTTT). where we use only the to- 
tal redshift distribution and do not split galaxies into red- 

9 Here, we fit a power-law with slope a minimizing the separation to 
all posted or- weighted points in the Q. m - cr 8 plane, and compute the ID 
marginalized mean of cr 8 (Q m /0.3)' r within Q m e [0.275,0.325]. 



shift bins. We find that the 2D and 3D analyses yield consis- 
tent results with substantially overlapping lcr regions, as ex- 
pected. Yet, the constraints from the 2D analysis shift towards 
lower cr 8 (An/0.3) 62 = 0.68 ± 0.11. The difference is not sur- 
prising given that the strongest contribution to the lensing sig- 
nal in COSMOS comes from massive str uctures near z ~ 0.7 
dScoville et alJ l2007b; M assev et alj|2007bl) . boosting the signal 
for high redshift sources, but leading to a lower signal for galax- 
ies at low and intermediate redshifts (see right panel of Fig.[8j. 
The 3D lensing analysis can properly combine these measure- 
ments, also accounting for the larger impact of sampling vari- 
ance at low redshifts. In contrast, the 2D lensing analysis leads 
to a rather low (but still consistent) estimate for cr 8 , due to the 
large number of low and intermediate redshift galaxies with low 
shear signal. 

The tomographic analysis also reduces the degeneracy be- 
tween D. m and <x 8 by probing the redshift-dependent growth of 
structure and distance-redshift relation, which differ substan- 
tially for a concordance ACDM cosmology and e.g. an Einstein- 
de Sitter cosmology (O m = 1). We summarize our parameter 
estimates in Table|2] also for the other cosmological models con- 
sidered in the following subsections. 

We also test our selection criteria for the optimized data vec- 
tor (Sect. l5.3l l by analysing several deviations from it for a flat 
ACDM cosmology. We find negligible influence if the small- 
est angular scales 6 < 0'5 or LRGs are included, suggesting that 
the measurement is robust regarding the influence of small-scale 
modelling uncertainties and intrinsic alignments between galaxy 
shapes and their surrounding density field. Performing the anal- 
ysis using only the usually excluded auto-correlations of the rel- 
atively narrow redshift bins 1 to 5, we measure a slightly lower 
cr 8 (Q m /0.3) ' 52 = 0.70 + 0.13, which is still consistent given 
the substantially degraded statistical accuracy. If intrinsic align- 
ments between physically associated galaxies contaminate the 
lensing measurement, we expect these auto-correlations to be 
most stro ngly affected. Howeve r, models predict an excess sig- 
nal (e.g. iHevmans et aLll2006bl) . whereas we measure a slight 
decrease within the statistical errors. Hence, we detect no sig- 
nificant indication for contamination by intrinsic galaxy align- 
ments. 

6.3.2. General (non-flat) ACDM cosmology 

We plot our constraints for a general ACDM cosmology without 
the assumption of flatness in Fig. [12] From the lensing data we 
find 

Q. A > 0.32 (90% conf.), 
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Fig. 12. Constraints on Q m , Qa, and <x 8 from our 3D weak lensing analysis of COSMOS for a general (non-flat) ACDM cosmology 
using our default priors. The contours indicate the 68.3% and 95.4% credibility regions, where we have marginalized over the 
parameters which are not shown. The non-linear blue-scale indicates the highest density region of the posterior. 



where our prior excludes negative densities Qa < 0. Based on 
our Q m - Qa constraints, we compute the posterior PDF for the 
deceleration parameter 

qo = -iia/a 2 = Q m /2 - Qa (22) 
as shown in Fig. [13] which yields 
qo < (96.0% conf.). 

Relaxing our priors to h = 0.72 + 0.08 (HST Key Project, 
iFreedman et al]l200ll). Q.uh 2 = 0.021 + 0.001 (Big-Bang nucle- 
osvnthesis. Ilocco et alJl2009h . and n s e [0.7, 1.2], weakens this 
constraint only slightly to 

qo < (94.3% conf., weak priors). 

Employing the recent distanc e ladder estimate 
h = 0.742 + 0.036 dRiess et alj [20091) instead of the HST 
Key Project constraint, we obtain qo < at 94.8% confidence. 

Our analysis provides evidence for the accelerated expansion 
of the Universe (qo < 0) from weak gravitational lensing. While 
the statistical accuracy is still relatively weak due to the limited 
size of the COSMOS field, this evidence is independent of ex- 
ternal constraints on Q m and Qa- 

We note that the lensing data alone cannot formally exclude 
a non-flat OCDM cosmology. However, the cosmological pa- 
rameters inferred for such a model would be inconsistent with 
various other cosmological probeo We therefore perform our 
analysis in the context of the well-established ACDM model, 
where the lensing data provide additional evidence for cosmic 
acceleration. 

6.3.3. Flat wCDM cosmology 

For a flat wCDM cosmology we plot our constraints on the (con- 
stant) dark energy equation of state parameter w in Fig.[14] show- 

10 For a lensing-only OCDM analysis the posterior peaks at 
Q. m - 0.1, erg = 1.4 (close to the prior boundaries). In the comparison 
with a ACDM analysis, the additional parameter Sl A c auses a penalty 
in the Bayesian model comparison (computed as in Kilbinger et al. 
2009b). This leads to an only slightly larger evidence for the non-fiat 
ACDM model compared to the OCDM model, with an inconclusive ev- 
idence ratio of 65:35. The evidence ratio becomes a "weak preference" 
(77:23) if we employ a (still conservative) prior cr 8 < 1. Hence, with 
this prior the ACDM model makes the data more than 3 times more 
probable than the OCDM model. 




qo 

Fig. 13. Posterior PDF for the deceleration parameter qo as com- 
puted from our constraints on Q m and Qa for a general (non-flat) 
ACDM cosmology, using our default priors (solid curve), and 
using weaker priors from the HST Key Project and Big-Bang nu- 
cleosynthesis (dashed curve). The line at qo = separates accel- 
erating (qo < 0) and decelerating (qo > 0) cosmologies. We find 
qo < at 96.0% confidence using our default priors, or 94.3% 
confidence for the weaker priors. 



ing that the measurement is consistent with ACDM (w = -1). 
From the posterior PDF we compute 

w < -0.41 (90% conf.) 

for the chosen prior w e [-2,0]. The exact value of this upper 
limit depends on the lower bound of the prior PDF given the non- 
closed credibility regions. We have chosen this prior as more 
negative w would require a worrisome extrapolation for the non- 
linear power spectrum corrections (Sect. 16.2] ). For comparison, 
we repeat the analysis with a much wider prior w e [-3.5,0.5] 
leading to a stronger upper limit w < -0.78 (90% conf.). While 
the COSMOS data are capable to exclude very large values 
w » -1, larger lensing data-sets will be required to obtain re- 
ally competitive constraints on w. 

To test the consistency of the data with ACDM, we com- 
pare the Bayesian evidence for the flat ACDM and wCDM 
m odels, which w e com pute in the PMC analysis as detailed 
in Kilbing er et alJ d2009bl) . Here we find completely inconclu- 
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Fig. 14. Constraints on O m and w from our 3D weak lensing 
analysis of COSMOS for a flat wCDM cosmology, assuming a 
prior w 6 [-2,0]. The contours indicate the 68.3% and 95.4% 
credibility regions, where we have marginalized over the param- 
eters which are not shown. The non-linear blue-scale indicates 
the highest density region of the posterior. 



Note that we did not apply this correction for the values given in 
the previous section and listed in Table |2 as we can only test it 
for the case of a flat ACDM cosmology. Additionally, we want 
to keep the results comparable to previous weak lensing studies, 
which we expect to be similarly affected. 

Having eliminated this last source of systematic uncertainty, 
we n ow estimate joint c onstraints with WMAP-5 CMB-only 
data (Dunklev et alJ 120091). c onducted similarly to the analy- 
sis by Kilbingereta D <|2009a). Here we assume a flat ACDM 
cosmology, completely relax our priors to Qb € [0.01,0.1], 
n s G [0.7, 1.2], h G [0.2, 1.4], and scale o"s for the lensing model 
calculation according to the Millennium Simulation results. 
Here we also marginalize over an additional 2% uncertainty in 
the lensing erg calibration to account for the dropped remaining 
mean shear calibration bias (0.8%, Sect. [3]) and limited accuracy 
of the employed residual shear correction (Sect.©, which we es- 
timate to be 1% in erg. From the joint analysis with WMAP-5 we 
find 



n - 266 +0 025+0 057 

li m - u.z,OO_ 023-0 .042 

+0.028+0.055 
0.029-0.060 



erg = 0.802 



(68.3%/95.4%conf., MS-calib.), 



sive probability ratios for wCDM versus ACDM of 52 : 48 
(w 6 [-2,0]) and 45 : 55 (w G [-3.5,0.5]), confirming that the 
data are fully consistent with ACDM. 

6.4. Model recalibration with the Millennium Simulation and 
joint constraints with WMAP-5 

Heitmann eFaXl d2008l) and iHilbert et ail (120091) found that the 
Smith et akl (120031) fitting functions slightly underestimate non- 



linear corrections to the power spectrum. To test whether this 
has a significant influence on our results, we performed a 
3D cosmological parameter estimation using the mean data 
vector of the 288 COSMOS-like ray-tracing realisations from 
the Millennium Simulation. Here we modify the strong pri- 
ors given in Sect. 16.11 to match the input values of the simula- 
tion (Q m = 0.25, o- 8 = 0.9, n s = 1, h = 0.73, Q b = 0.045), and 
find cr 8 = 0.947 ± O.OOffi for Q m = 0.25. This confirms the re- 
sult of iHeitmann etal] 120081) and IHilbert et al.1 (l2009t) . indi- 
cating that models based on ISmifh et all (120031) slightly under- 
estimate the shear signal, hence a larger o"8 is required to fit 
the data. Here we use actual reduced shear estimates from the 
simulation, but employ shear predictions, as done for the real 
data (see Sect.|4]l. Using shear estimates from the simulation 
yields o"8 = 0.936 + 0.006. Hence, a minor contribution to the 
overestimation of o- 8 is c aused by the negligence of reduced 
shear corrections (see also iDodelson et alj|2006t IShapiroi r2009; 
iKrause & Hiratall2009h . 

To compensate for this underestimation of the model pre- 
dictions and reduced shear effects, we scale our derived con- 
straints on o"8 for a flat ACDM cosmology by a factor 
0.9/0.947 a 0.95(0, yielding 



o- 8 (Q m /0.3)' 



0.51 



0.75 + 0.08 (68.3% conf ., MS-calib.). 



1 1 Here we have scaled the uncertainty for the mean ray-tracing data 
vector from the uncertainty for a single COSMOS-like field assuming 
that all realizations are completely independent. This is slightly opti- 
mistic given the large but finite volume of the simulation, and fact that 
the realizations were cut from larger fields. 

12 We expect that this correction factor depends on cosmological pa- 
rameters. Yet, considering the weak lensing degeneracy for Q m and cr 8 , 
the input values of the Millennium Simulation are quasi equivalent to 



which reduces the size of WMAP-only lo" (2o-) error-bars on 
average by 21% (27%). We plot the joint and individual con- 
straints in Fig. [15] illustrating the perfect agreement of the two 
independent cosmological probes. 



CO 

b 



1.4 



1.2 



1.0 



0.8 



0.6 



0.4 



0. 



s 



1 / \\ 



Lensing+WMAP 
Lensing 
WMAP 




0.2 



0.4 



0.6 



0.8 



1.0 



m 



Fig. 15. Comparison of the constraints on Q. m and o"8 for a 
flat ACDM cosmology obtained with our COSMOS analysis 
(dashed), WMAP-5 CMB data (dotted), and joint constraints 
(solid). The contours indicate the 68.3%, 95.4%, and 99.7% 
credibility regions. Note that the weak lensing alone analysis 
uses stronger priors. The weak lensing constraints on o"8 have 
been rescaled to account for modelling bias of the non-linear 
power spectrum and reduced shear corrections according to the 
ray-tracing constraints from the Millennium Simulation. 



cr 8 ^ 0.82 for fi m = 0.3, which is sufficiently close to our constraints to 
justify the application. 
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7. Summary, discussion, and conclusions 

We have measured weak lensing galaxy shear estimates from the 
HST/COSMOS data by applying a new model for the spatially 
and temporally varying ACS PSF, which is based on a principal 
component analysis of PSF variations in dense stellar fields. We 
find that most of the PSF changes can be described with a single 
parameter related to the HST focus position. Yet, we also correct 
for additional PSF variations, which are coherent for neighbour- 
ing COSMOS tiles taken closely in time. We employ updated 
parametric corrections for charge-transfer inefficiency, for both 
galaxies and stars, removing earlier modelling uncertainties due 
to confused PSF- and CTI-induced stellar ellipticity. Finally, we 
employ a simple correction for signal-to-noise dependent shear 
calibration bias, which we derive from the STEP2 simulations of 
ground-based weak lensing data. Tests on simulated space-based 
data confirm a relative shear calibration uncertainty \m\ < 2% 
over the entire used magnitude range if this correction is applied. 
We decompose the measured shear signal into curl-free E-modes 
and curl-component B-modes. As expected from pure lensing, 
the B-mode signal is consistent with zero for all second-order 
shear statistics, providing an important confirmation for the suc- 
cess of our correction schemes for instrumental systematics. 

We combine our shear catalog ue with excelle nt ground- 
based photometric redshifts from lllbert et ail ((2009) and care- 
fully estimate the redshift distribution for faint ACS galaxies 
without individual photo-zs. This allows us to study weak lens- 
ing cross-correlations in detail between six redshift bins, demon- 
strating that the signal indeed scales as expected from General 
Relativity for a concordance ACDM cosmology. 

We employ a robust covariance matrix from 288 simu- 
lated COSMOS-like fields o btained from ray-tra cing through 
the Millennium Simulation dHilbert et al.l 120091) . Using our 
3D weak lensing analysis of COSMOS, we derive con- 
straints cr 8 (Q m /0.3) - 51 = 0.79 + 0.09 for a flat ACDM cos- 
mology, us in g non -linear power spectrum corrections from 
Smit h et al.l Q2003). A recalibration of these predictions 
based on the ray-tracing analysis changes our constraints to 
cr 8 (Q m /0.3) a51 = 0.75 + 0.08 (all 68.3% conf.). Our results are 
perfectly consistent with WMAP-5, yielding joint constraints 
f2 ra = 0.266!°;°!!°;^, <r g = 0.802+°™;°*; (6g 3% and g5 4% 
confidence). They also agree with weak lensing results from 
the CFHTLS-W ide dFu et al 1 120081) and recent galaxy cluster 
constraints from iMantz e t al. (2009) within lcr. Our errors in- 
clude the full statistical uncertainty including the non-Gaussian 
sampling variance, Gaussian photo-z scatter, and marginaliza- 
tion over remaining parameter uncertainties, including the red- 
shift calibration for the faint i + > 25 galaxies. 

Our results are consistent with the 3D lensing con- 
straints cr^tlJO.l)^ = 0.866 + 0.033 (stat.)±°;°§ ( sys .) from 
Mass ev et al.l (l2007d) assuming non-lin ear power spectrum cor- 
rections according to Smith et al. (2003), at the ~ lcr level. The 
analyses differ systematically in the treatment of PSF- and CTI- 
effects, where the success of ou r methods is confirm ed by the 
vanishing B-mode. Furthermore, Mas sey et al.1 (l2007d) emplo y 
earlier photo-zs based on fewer bands (iMobasher et al.1 12007). 
Note that the analysis of Ma ssev et al.l (l2007d) yields tighter sta- 
tistical errors, which may be a result of their covariance estimate 
from the variation between the four COSMOS quadrants. This 
potentially introduces a bias in th e covariance inver sion due to 
too few independent realisations (lHartlap et alJ2007h . While the 
absolute calibration accuracy of the shear measurement method 
was estimated to be the dominant source of uncertainty in their 
error budget, we were able to reduce it well below the statis- 



tical error level. As a further difference, our analysis employs 
photometric redshift information to reduce potential contamina- 
tion by intrinsic galaxy alignments, where we exclude the shear- 
shear auto-correlations for the relatively narrow redshift bins 
1 to 5 to minimize the impact of physically associated galax- 
ies. In addition, we exclude luminous red galaxies, which were 
found to carry the strongest intrinsic alignment with the density 
field of their large- scale structure environment causing the shear 
(Hira ta et al.l 120071) . Finally, we do not include angular scales 
9 < 0'5 due to increased modelling uncertainties for the non- 
linear power spe ctrum. 

Similarly to Masse v et al.l d2007d) . we obtain a lower esti- 
mate a- 8 (Q m /0.3) 62 = 0.68 ± 011 for a non-tomographic (2D) 
analysis, assuming Smi th et all (|2003) power spectrum correc- 
tions. The lower signal compared to the 3D lensing analysis is 
expected, given that the mos t massive structures i n COSMOS 
are located at 0.7 < z < 0.9 (IScoville et al.l l2007bl) . creating a 
strong shear signal for high redshift sources only, which is de- 
tected by the 3D analysis. In contrast, the bulk of the galaxies 
in the 2D lensing analysis are located at too low redshifts to be 
substantially lensed by these structures, yielding a relatively low 
estimate for cr 8 . Nonetheless, as sampling variance is properly 
accounted for in our error analysis, the constraints are still con- 
sistent. 

For a general (non-flat) ACDM cosmology, we find a neg- 
ative deceleration parameter qo < at 96.0% confidence us- 
ing our default priors, and at 94.3% confidence if only priors 
from the HST Key Project and BBN are applied. Hence, our 
tomographic weak lensing measurement provides independent 
evidence for the accelerated expansion of the Universe. For a 
flat wCDM cosmology we constrain the (constant) dark energy 
equation of state parameter to w < -0.41 (90% conf.) for a prior 
w e [-2, 0], fully consistent with ACDM. Our dark energy con- 
straints are still wea k compared to recen t results from inde- 
pendent probes (e.g. iKowalski et al.l 12008b iHicken et aD |2009: 
[ Allen et alj200llMantz et al.ll2008l 120091: 1 Vikhlinin e7ail2 009: 
iKomatsu et al.ll2009l) . This is solely due to the limited area of 
COSMOS, leading to a dominant contribution to the error bud- 
get from sampling variance. 

While the area covered by COSMOS is still small 
(1.64deg 2 ), the high resolution and depth of the HST data al- 
lowed us to obtain cosmological constraints which are compa- 
rable to results from substantially larger ground-based surveys. 
However, note that HST was by no means designed for cos- 
mic shear measurements. In contrast, future space-based lensing 
mission such as EuclicQ or JDEM^ will be highly optimised 
for weak lensing measurements. High PSF stability, a much 
larger field-of-view providing thousands of stars for PSF mea- 
surements, carefully designed CCDs which minimize charge- 
transfer inefficiency, and improved algorithms will remove the 
need for some of the empirical calibrations employed in this pa- 
per. 

In order to fully exploit the information encoded in the weak 
lensing shear field, second-order shear statistics, as used here, 
can be complemented with higher-order shear statisti cs to probe 
the non-Gaussianity of t he matter distribution (e.g. Berge et al. 
120091: |Vafaeietal.ll2010l). Ba sed on our COSMOS shear cata- 
logue, |Se^Solonretar| {2010) present such a cosmological anal- 
ysis using combined second and third-order shear statistics. 

Finally, we stress that weak lensing can only provide preci- 
sion constraints on cosmological parameters if sufficiently ac- 

13 http://sci.esa.int/euclid 

14 http://jdem.gsfc.nasa.gov/ 
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curate models exist to compare the measurements to. Our anal- 
ysis of the relatively small COSMOS Survey is still limited by 
the statistical measurement uncertainty, for which our approx- 
imate model recalibration using the Millennium Simulation is 
sufficient. Most of the cosmological sensitivity in COSMOS 
comes from quasi-linear and non-linear scales. We cut our anal- 
ysis only at highly non-linear scales 9 < 0'5, corresponding to a 
comoving separation of ~ 360 kpc at z = 0.7 (roughly the red- 
shift of the most massive structures in COSMOS). At these 
scales non-linear power spectrum corrections have substantial 
uncertainties, in particular due to the influence of baryons (e.g. 
iRudd et aU2008b . Given that our results are basically unchanged 
if even smaller scales are included (insignificant increase in cr 8 
by < 1%), we expect that the model uncertainty for the larger 
scales should still be sub-dominant compared to our statistical 
errors. However, analyses of large future surveys will urgently 
require improved model predictions including corrections for 
baryonic effects, also for dark energy cosmologies with w + — 1, 
and optionally also for theories of modified gravity. Once these 
are available, careful analyses of large current and future weak 
lensing surveys will deliver precision constraints on cosmologi- 
cal parameters and dark energy properties. 
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Appendix A: Additional image calibrations 

In this appendix we describe additional calibrations which 
we apply to the flat-fielded _/?f images before running 
MultiDrizzle. 

Background subtraction. We perform a quadrant-based back- 
ground subtraction due to an anomalous bias level variation be- 
tween the four ACS read-out a mplifiers. Here we dete ct and 
mask objects with SExtractor (Berti n & Arnoutslll99 6l. com- 
bine this mask with the static bad pixel mask, and estimate 
the background as the median of all non-masked pixels in the 



Table A.l. Lower and upper cr thresholds for pixel masking with 
ccdmask in the bias variance and the sky-subtracted and object- 
masked median images. 



Image type 


lsigma 


hsigma 


Bias variance 


100 


25, +5 if more than 2.5% masked 


Median, gain=l 


13 


11 


Median, gain=2* 


15 


15 



*: The COSMOS images were taken with gain = 1, whereas the 
gain = 2 setting has been applied for some of the HAGGLeS fields 
(Marshall et al. 2010, in prep.). While we do not include these fields 
in the current analysis, they have been processed with the same pipeline 
upgrades described here. Hence, we list these values for completeness. 



quadrant. We modulate the offset from the mean background 
level with the normalised inverse flat-field to correct for the fact 
that the improperly bias-subtracted image has already been flat- 
fieldedB 

Bad pixel masking. We manually mask satellite trails and scat- 
tered stellar light if its apparent sky position changes between 
different dither positions, allowing us to recover otherwise unus- 
able sky area. In addition, we update the static bad pixel mask 
rejecting pixels if: 

- their dark current exceeds 0.04 e~/sec in the associated dark 
reference file (default 0.08 e~/sec), or 

- they are affected by variable bias structures, which we iden- 
tify in a variance image of five subsequent bias reference 
frames taken temporally close to the science frame consid- 
ered, or 

- they show significantly positive or negative values in a me- 
dian image computed from 50 background-subtracted and 
object-masked COSMOS frames taken closely in time, in- 
dicating any other semi-persistent blemish. 

The latter two masks mainly aim at the rejection of variable bias 
structures which show up as positive or negative bad column 
segments in the stacked image if not properly masked. For the 
mask creation we utilize the IRAF task noao . imred . ccdmask. 
It computes the local median signal and rms variation in moving 
rectangles. A pixel is then masked if its values is either lsigma 
below or hsigma above the local median value. This is done for 
individual pixels and sums of pixels in column sections, where in 
the latter case the background dispersion is scaled by the square 
root of the number of pixels in the section. Finally each col- 
umn is scanned for short segments of un-flagged pixels in be- 
tween masked pixels. We additionally mask these segments if 
their length is less than 15 pixels. We summarise the values ap- 
plied for the thresholds lsigma and hsigma in Table |A~T1 Due 
to variations in image noise properties they do not perform opti- 
mally in all cases, so that we iteratively increase hsigma by +5 
if otherwise more than 2.5% of the pixels in the bias variance 
image would be masked. 



15 This procedure performs well for relatively empty fields such as 
the large majority of the COSMOS tiles. For fields dominated by a very 
bright star or galaxy, it can, however, lead to erroneous jumps in the 
background level. Hence, we generally adopt a maximal accepted dif- 
ference in the background estimates of 4e~, which, if exceeded, leads 
to a subtraction of the minimum background estimate for all quadrants. 
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Noise model. We compute a rms noise model for each pixel as B.2. Tests with simulated space-based data 



ERR = F~ 



4 



sF + f ex pl> + erf 



+ 7 2 V [e"], 



(A.l) 



with the normalised flat-field F, the sky background s [e ], the 
dark reference frame D [e~/s], the exposure time f exp [s], the 
read-noise cr r 5 e~, and the bias variance image V [counts 2 ] 
described in the previous paragraph, which requires scaling with 
the gain y [e~ /count]. Containing all noise sources except object 
photon noise, this rms model is used for optimal pixel weighting 
in MultiDrizzle. 



Appendix B: Correction for PSF and CTI effects 

B.1. Summary of our KSB+ implementation 

We measure galaxy shapes using the Erben et al. I (l200Th im- 
plementation of the KSB+ formalism (jKaiser et akl Il995t 
Lup pino & Kaiseil 1 1997t iHoekstraet ail 1 1998b. as done in the 
earlier ACS weak lensing analysis of Schrabba ck et al.l (120071) . 
Object ellipticitie{3 



e = e\ + 1^2 



fin -622 + 21612 



(B.l) 



<2n + 622 

are measured from weighted second-order brightness moments 

Q u = j d 2 6 W rt mWm, hje{h2), (B.2) 

where W r is a 2D Gaussian with dispersion r g . The response 
of a galaxy ellipticity to reduced gravitational shear g and PSF 
effects is given by 



(B.3) 



with the (seeing convolved) intrinsic source ellipticity e s and the 
"pre-seeing" shear polarisability 



a/5 



a/5 



[(P 



yd 



U p 



(B.4) 



where the shear and smear polarisability tensors P sh and P sm 
are calcu lated from hig h er-ord er brightness moments as de- 
tailed in iHoekstra et alJ dl998l) . The PSF anisotropy kernel 
q* a = (P sm *)~pe* and ratio of P sh * and P sm * must be measured 
from stars and interpolated for each galaxy position, where we 
approximate the latter as T* = Tr [f sh *] /Tr [P sm *]. 

In the application of the KSB+ formalism several choices 
lead to su btle differences between different KSB implementa- 
tions, see iHevmans et alJ (l2006al) for a detailed comparison. 
In short, we use sub-pixel interpolation for integral evalua- 
tions, measure galaxy shapes with r g = r f , the SExtractor flux- 
radius, and apply PSF measurements computed with the same 
filter scale as used for the corresponding galaxy (interpolated 
between 24 values with 1 < r„ < 15 pixels). We invert the P g 
tensor as measured from individual galaxies using the approx- 
imation (P g )~ l = 2/Tr\P g ] comm only applied to reduce noise 
dErben et alj|2001l) . In contrast to Schrabb ack et alJ d2007l) we 
do not apply a constant calibration correction, but employ the 
signal-to-noise dependent correction (0. 



16 We adopt the widely used term "ellipticity" here, but note that, 
strictly speaking, l lB. It corresponds to the definition of the polarisation. 



We test our KSB+ shape measurement pipeline on simu- 
lated space-based weak lensing data with ACS-like proper- 
ties, which were provided for testing in the framework of the 
Sh ear Testing Program mJ^. The images were created with 
the iMassey et al.l d2004j) image simulation pipeline, which uses 
shapelets (Refregie r & Baconl 120031: IMassev & Refregierl l2005) 
to model galaxy and PSF shapes, as alrea dy employed for 
the STEP2 simulations dMassev et al.ll2007al) . All images have 
4k x 4k pixels of size 0'.'04, HST-like resolution, and a depth 
equivalent to 2ks of ACS imaging. The data are subdivided into 
eight sets with different PSFs (|e*| < 7%), seven of which uti- 
lize TinyTirr[3 ACS PSF models, and one was created by stack- 
ing stars of similar ellipticity in an ACS stellar field (M). One 
of the sets uses simplified exponential profiles for galaxy mod- 
elling (F), while the others include complex galaxy morpholo- 
gies modelled with shapelets. Four sets comprise 100 images, 
while the others include 200 frames. Within each set, the images 
are split into "rotated pairs", where the intrinsic galaxy elliptic - 
ities in one frame resemble 90 d egree-rotated version s from the 
other frame, an approach used in Mass ev et al.l (l2007ab to reduce 
the analysis uncertainty due to shape-noise. Galaxies are sheared 
with \g\ < 0.06 and convolved with the PSF, both effects being 
constant within one frame, but with varying g within one set. 
Realistic image noise was added similarly to the STEP2 analy- 
sis, except that no noise correlations were introduced. 

We analyse the images with the same pipeline and cuts as 
the real COSMOS data, with the only difference that the PSF is 
assumed to be constant across the field, but still measured from 
the simulated stars. Fig. IB. II shows the mean calibration bias m 
and PSF anisotropy residuals c defined in (0, separately for each 
image simulation set, estimated from match ed galaxy pairs (for 
details on this fit see IMassev eta l. 2007a). While some data- 
sets deviate from the optimal m — c — 0, the residuals are at a 
level which is negligible compared to the statistical uncertainty 
of COSMOS. Combining all sets and both shear components, 
we estimate the mean calibration bias m = +0.008 + 0.002, and 
the scatter of the PSF anisotropy residuals cr c = 0.0006. 

As discussed in Sect. [3] a possible magnitude-dependence of 
the shear calibration bias m is particularly problematic for 3D 
weak lensing studies. We therefore study m as a function of mag- 
nitude in Fig. IB. 21 both for the simulated ground-based STEP2 
and the simulated space-based ACS-like data. Although the ap- 
plied correction (|6]l was determined from the simulated ground- 
based data, it also performs very well for the ACS-like simula- 
tions, showing its robustness. Over the entire magnitude range 
the remaining calibration bias is \m\ < 0.02, which is negligible 
compared to our statistical errors. 



B.3. Stellar fields 

We have analysed 700 i%u exposures of dense stellar fields, 
which were taken between 2002 Apr 18 and 2006 Jun 03 and 
contain at least 300 non-saturated stars with S /N > 50 (for 
r g = 1.5 pixels). This large set enables us to study in detail 
the impact of CTI on stars, as well as the temporal and posi- 
tional ACS PSF variation, which cannot be achieved from the 
COSMOS exposures due to their low stellar density. 

We determine both CTI and PSF models for the cosmic ray- 
cleansed COR images before resampling, and their resampled 



http : //www . physics . ubc . ca/~ heymans/step . html 
http : //www . stsci . edu/sof tware/tinytim/ 
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Fig. B.l. Shear calibration bias m and PSF anisotropy residuals c 
as measured in the simulated ACS-like lensing data. The left and 
right panels show the results for the j\ and j2 shear components 
respectively. Each letter corresponds to a different PSF model. 
Although some data-sets deviate from the optimal m = c = 0, 
the residuals are at a level which is negligible compared to the 
statistical errors for COSMOS. 
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Fig. B.2. Magnitude-dependence of the shear calibration bias 
for our KSB implementation after correction for S /W-dependent 
bias according to ©. The top panel shows r esults for the STEP 2 
simulations of ground-based lensing data (Massev et al. 2007a), 
which have been used to derive ©, where we have excluded 
the untypically elliptical PSFs D and E. The bottom panel shows 
the remaining calibration bias for the ACS-like simulations of 
space-based lensing data. In both panels we plot the average 
computed from all PSF models and the two shear components, 
with error-bars indicating the uncertainty of the mean. Despite 
the very different characteristics of the two sets of simulations, 
(O performs also very well for the ACS-like data, with a bias 
\m\ < 0.02 over the entire magnitude range. The remaining cali- 
bration uncertainty is negligible compared to the statistical errors 
for COSMOS. 



(but not stacked) counterparts (DRZ). The reason is that resam- 
pling unavoidably adds extra noise, hence it is best to fit the 
available stars in a galaxy field exposure before resampling. Yet, 
the combined PSF model for a stack has to be determined from 
resampled image models according to the relative dithering. For 
the CCW-image analysis we employ a fixed Gaussian filter scale 
r g = 1.5 pixels, in order to m aximize the fitting signal-to-noise 
(see Schrabback et al. 2007), and characterize the PSF by the 
ellipticity <?* and stellar half-light radius r^ as suggested by 



iJee et al For the DRZ images we require CTI-corrected 

PSF models for all 24 values of r„ used for the galaxy correction. 

B.4. Stellar CTI correction 

CTI charge trails stretch objects in the readouts-direction, lead- 
ing to an addit ional negative e\ el lipticity component. Internal 
calibrations ( Mutchj er & SManiiil 120051) . photometric studies 



(e.g. IChiaberge et alj 120091) . as well as the a nalysis of warm 
pixels (Masse v et al.l 120101) and cosmic rays (IJee et al.l [2009) 



demonstrate that the influence of CTI increases linearly with 
time and the number of y-transfers, where the latter has also been 
shown for the influence on galaxy ellipticities by R hodes et al.l 
(2007). In addition, the limited depth of charge traps leads to a 
stronger influence of CTI for faint sources, which lose a larger 
fraction of their charge than bright sources. Likewise, the effect 
is reduced for higher sky background values leading to a fraction 
of continuously filled traps. Here we only study the effect of CTI 
on stars, whereas galaxies will be considered in App. IB.6l 

Following Chiaberge et a l. (2009), we assume a power-law 
dependence on sky background and integrated flux as measured 
in apertures of 4.5 rf 5.8 pixels, leading to the parametric CTI 
model 



cti,*/ \ f 

e l (r s ) = -e 



/ FLUX \ 8 I SKY \ s l t \ 
l(rg Hl0 4 e- j Ue ] llOOOdJ 

1 Jtrans \ 



V2048/ 



(B.5) 



with the time t = MJD - 52340 since the installation of ACS on 
2002 Mar 08, and the number of y-transfers y tmns . We expect 
that the normalisation e ( [(r g ) and power law exponents F(r g ) and 
S(r g ) depend on the Gaussian filter scale r g of the KSB ellip- 
ticity measurement. E.g., for a measurement of the PSF core 
with small r g , charge traps may already be filled by electrons 
from the outer stellar profile, leading to an expected strong flux- 
dependence. On the contrary, the PSF wings measured at large r„ 
will be more susceptible to trap filling by background electrons, 
leading to a stronger sky-dependence. 

In order to separate CTI and PSF effects we make use of 
the fact that CTI-induced ellipticity is expected to depend on 
flux, while PSF ellipticity is flux-independent. In our analy- 
sis of stellar field exposures we first fit the spatial ellipticity 
variation of bright non-saturated stars with S /N > 50 using a 
third-order polynomial in each chip, and apply this model to 
all stars with S/N > 5. For the high S/N stars used in the fit, 
the strongest ellipticity contribution comes from the spatially 
varying PSF. Yet, for these stars the polynomial fit also corrects 
for the position-dependent but flux-averaged CTI effect, lead- 
ing to a net negative e\ ellipticity for fainter than average stars 
(CTI under-corrected), and net positive e\ for brighter stars (CTI 
over-corrected). For even fainter stars with S/N < 50 we expect 
an increasingly more negative e\ ellipticity component. Hence, 
the CTI influence can be measured from the flux-dependence 
of the polynomial-coiTected residual ellipticity, as illustrated in 
Fig. lB.3l for four example exposures. We note a turnaround in the 
CTI flux-dependence for some exposures at low S/N ~ 5 - 10 
(right panels in Fig. IB. 31 ), which was also reported for CTI mea- 
surem ents from cosmic rays and further investigated bv lJee et ail 
d2009l) . This does not affect our stellar models, given that we 
only use SJN > 20 stars both for PSF measurement and to con- 
strain ( IB.5l l. Yet, it suggests that CTI models may not be valid 
over very wide ranges in signal-to-noise, motivating the use 
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Fig. B.3. CTI-induced stellar ellipticity for four example stel- 
lar field exposures: The bold points show the mean stellar e\ 
ellipticity-component as a function of stellar flux after subtrac- 
tion of a spatial third-order polynomial model derived from 
bright stars (S /N > 50) to separate PSF and CTI effects. Each 
stellar ellipticity has been scaled to a reference number of y tmn s = 
2048 parallel readout transfers. The curves show the CTI model 
( IB. 51 , where the fit parameters have been jointly determined 
from S IN > 20 stars in all 700 exposures, and an offset shown by 
the horizontal line has been applied, corresponding to the mean 
CTI model ellipticity of the bright stars used for the polynomial 
interpolation. The crosses indicate the corrected ellipticities af- 
ter subtraction of the CTI model. Note the strong increase of 
the CTI-induced ellipticity with time (top left to bottom right) 
and moderate dependence on sky background {top right versus 
bottom left at similar times). Also note the turnaround occurring 
for faint stars at ~ 1 000e~ (co rresponding to S/N ~ 5 - 10) in 
the right panels, see iJee et al .1 (120091) for a further investigation 
of this effect. The plots shown correspond to the non-resampled 
CCW-images with ellipticities measured using a Gaussian filter 
scale of r„ — 1.5 pixels. 



of a separate model for the typically much fainter galaxies in 
App.ESl 

We determine the three fit parameters in ( 1B.51 I jointly from 
the polynomial-corrected residual ellipticities in all stellar expo- 
sures. For each exposure it is necessary to add an offset, which 
has been linearly scaled with y tmns for each star, in order to com- 
pensate for the flux-averaged correction included in the poly- 
nomial fit. We compute this offset within the non-linear fitting 
routing f° r a given set of fit parameters from the positions and 
fluxes of the bright stars used in the polynomial fit, and apply it 
to all stars. 

We conduct this fit both for the COR-image 
ellipticities (r g = 1.5) yielding best fitting values 
(e° F,S) = (0.0073 ± 0.0002, 0.65 ± 0.02, 0.06 ± 0.01), and 



for the resampled DRZ-images for all values of r g . For the latter 



19 For the non-linear CTI fits we utilize the CERN Program Library 
MINUIT: http://wwwasdoc .web.cern. ch/wwwasdoc/minuit/ 
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Fig. B.4. Dependence of the best fitting parameters of the stellar 
CTI model ( IB. 51 ) on the Gaussian filter scale r g used for shape 
measurements in the resampled DRZ-images. The curves corre- 
spond to the fitting functions (1B.6I I . 

Table B.l. Fitted coefficients for the r g -dependent CTI- 
ellipticity model ( IB.6l l in the resampled DRZ frames. 
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we adjust the S /N cuts in order to keep enough stars for large 
r g . The best fit values are shown in Fig. IB.4l as function of r g , 
indeed confirming the expected trends. We provide the fitting 
functions 



to 



;=3 

Z 

7=0 



F(r g ) 



7=3 
7=0 



S(r g ) 



7=3 
7=0 



(B.6) 



where the coefficients are listed in Table lB.ll being valid for 
1 < r g < 15 pixels. We correct the ellipticities of all stars both 
in the stellar and galaxy fields with the derived models, as im- 
plicitly assumed in the following sections. 



B.5. Principal component correction for the time-dependent 
ACS PSF 

As discussed in Sect.[3] ACS PSF variations are expe cted to be 
mostly caused by changes in telescope f ocus (e.g. iKrisl 12003: 
lLallo et aLll2006t [Anderson & Kingll2006h . If the temporal vari- 
ations indeed depend on one physical parameter only, it should 
be possible to construct a one-parametric PSF model, which 
can be well constrained with the ~ 10-20 stars available in an 
ACS field at h igh galactic l a titude s. Such an approach was im- 
plemented bv iRhodes et al.l ([2007 ), who measure the mean fo- 
cus offset for a COSMOS stack from simulated focus-dependent 
TinyTim PSF models. They then interpolate the ACS PSF be- 
tween all stars in COSMOS using polyn omial functions depen- 
dent on both positi on and focus offset (Leauth aud et al.l [2007: 
iMassev et al.l l2007d) . However, as sugg ested by the residua l 
aperture mass B-mode signal found by Masse vet al.l d2007cl) . 



Tim Schrabback et al.: Evidence for the accelerated expansion of the Universe from 3D weak lensing with COSMOS 



21 




5 10 15 

Principal component 1 



Fig. B.5. Fractional PCA eigenvalues for the PSF variation in 
700 z'gi4 ACS stellar field exposures. The dashed (solid) curve 
has been computed considering the variation of e\ and e* 2 (e\, e* 2 , 
and r*). The dominant first principal component contains 97% 
(95%) of the variation and is caused by focus variations. 

this approach appears to be insufficient for a complete removal of 
systematics. In an alternative approach, Schrabb ack et alj d2007l) 
fit the stars present in each galaxy field exposure using a large 
library of stellar field PSF models. While this approach led to no 
significant residual systematics within the statistical accuracy of 
GEMS, it is also not sufficient for the analysis of the much larger 
COSMOS data-set. We therefore implement a new PSF interpo- 
lation scheme based on principal component analysis. It effec- 
tively combines the idea of exposure-based empirical models, 
which optimall y account for time variations and relative dither- 
ing dSchrabbac k et al. I l200l . with the aim to describe the PSF 
var iation with a s i ngle p arameter dRhodes et al. 2007). 

iJarvis & Jain! (2004) introduced the application of principal 
component analysis (PCA) for ground-based PSF interpolation, 
which we adapt here to obtain well-constrai ned PSF mode ls for 
our ACS weak lensin g fields. Note that iJee et al.l d2007l) and 
iNakaiima et alj d2009l) employed PCA to efficiently describe the 
two-dimensional ACS PSF shape, which they then spatially in- 
terpolated with normal polynomial functions. Thi s is conceptu- 
ally ve ry different to the approach suggested by IJarvis & Jain! 
(2004) and used here, which employs PCA for the spatial and 
temporal interpolation of certain quantities needed for PSF cor- 
rection, such as the stellar ellipticity e*. 

We represent all quantities which we want to interpolate as 
p a . This includes e* v el,r^ measured in the COR images for 
r g = 1.5 pixels, but also e*, el, q\, q\, T* as measured in the DRZ 
images for varying r„. The only exception is when we specif- 
ically allude to COR quantities, which only includes the first 
group. 

The first step of the PCA analysis is to fit the positional vari- 
ation of the three COR PSF quantities in all j <N = 700 stellar 
exposures jointly for both chips with 3rd-order polynomials 

m 

yielding m = 10 coefficients each, where we generally denote 
polynomials using a capital P with the order indicated by the 
superscript. Here we account for the gap between the chips and 
rescale the pixel range to the interval x,y e [0, 1 ] . While this fit is 



unable to account for some small-scale features such as a small 
discontinuity between the two chips, it captures all major large- 
scale PSF variations and is very well constrained by the required 
> 300 stars. For each exposure we arrange the M = 3 X m — 30 
polynomial coefficients in a data vector dj, with components 
(now i < M). We then subtract the mean vector and divide each 
component with an adequately chosen normalization n,, yielding 
the modified data vector dj with 

da - dj 

dij = . (B.8) 

Hi 

We then aiTange all modified data vectors into a M x N dimen- 
sional data matrix D - \d\, The central step of the PCA 
is a singular value decomposition D = WLV 1 , where the or- 
thonormal matrix W consists of the singular vectors of D, and 
the diagonal matrix T, contains the ordered singular values su of 
D as diagonal elements. Here the /th largest singular value cor- 
responds to the /th singular vector, which is also named the /th 
principal component. In the coordinate system spanned by the 
singular vectors, the matrix C = DD T = WL1, T W T = WAW T 
corresponding to the covariance matrix for n, ■ - 1, becomes di- 
agonal, where the sorted eigenvalues A/ - sf, are equal to the 

variance of the vectors dj along the direction of the /th principal 
component. 

Note that the relative values and absolute scale of the eigen- 
values Ai depend on the normalisations n,. Uniform n, = 1 would 
not be adequate given that we combine PSF quantities with dif- 
ferent units (dimensionless e* versus rjj in pixels). A congela- 
tion analysis with n,- = <x ; = ~~ di) 2 /N could be used, 

but here relatively stable polynomial coefficients with small cr, 
would unnecessarily add noise, effectively increasing the rela- 
tive eigenvalues of higher principal components. Aiming at a 
compact description of most of the actual PSF variation in the 
field with a small number of important principal components, 
we employ the normalisation 

n i = Qi i + l)(y i + l)o- a , (B.9) 

where we use the mean variance of all coefficients belonging to 
the corresponding PSF quantity p a : 

i W j n _ 

Z fit**-*?- ( B - 10) 

"=<mi„,„ /=1 

In this way the c/, ; in dB.8t become dimensionless, all three PSF 
quantities contribute similarly to the total variation, and the un- 
desired noise from relatively stable polynomial coefficients is 
avoided, as their variation is averaged with that from the less 
stable coefficients. The pre-factor in ( IB.9b equals the inverse of 

the integral J Q di dy xry 1 of the corresponding polynomial term 
in ( IB.7b . It accounts for our aim to scale according to the actual 
PSF variation, where e.g. a Oth-order term affects the whole field 
while a 3rd-order term with similar amplitude gets lower weight 
as it contributes substantially in a smaller area only. 

We plot the fractional eigenvalues in Fig. IB. 51 once using the 
analysis as described above (solid curve) and once considering 
only the two ellipticity components without (dashed curve). In 
both cases the first principal component is clearly dominant, con- 
tributing with 95% (97%) of the total variance. We identify this 
variation as the influence of focus changes, which are expected 
to dominate the actual PSF variation. The reason why the second 
principal component has a larger eigenvalue if rj| is included in 
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Fig. B.6. Variation of 700 /8i4 stellar field exposures in the space spanned by the first three principal components, which have been 
computed using the polynomial coefficients of e\, el, and rjj. Note the different axis scales. The non-linear dependence in the left 
panel is caused by the different response of PSF ellipticity and size on defocus, and leads to the increased eigenvalue Xi in Fig. IB. 51 
if the variation is included in the PCA. The data points have been split according to the velocity aberration plate scale factor 
VAF ACTOR. The fact that the three subsets scatter differently for fixed y\j shows that deviations from pure focus variations are not 
completely random, but depend on orbital parameters and may hence be coherent for surveys such as COSMOS. 



the analysis (fractional A2 = 1 . 1 % versus A2 = 0.6%) can be seen 
if we project the data variation onto the space spanned by the 
singular vectors Y = W T D, with components yy. Looking at the 
yij - yij variation in the left panel of Fig. IB. 61 where has been 
included, we see that the data points roughly follow a quadratic 
curve in the plane defined by the first two singular vectors. The 
reason for this is the linear response of PSF ellipticity on defocus 
caused by astigmatism, w hile PSF width res ponds to leading or- 
der quadratically (see e.g. lJarvis et alj |2008). Given that PCA is 
a purely linear coordinate transformation, it is not capable to di- 
rectly capture this one-parametric variation (separation between 
primary and secondary mirror) with a single principal compo- 
nent. This is only possible if PSF quantities with the same de- 
pendence on physical parameters are included, hence the smaller 
A2 if only the two ellipticity components are considered. Thus, 
for other applications it might be more favourable to perform 
a PCA analy sis for each consider ed PSF quantity separately, as 
also done bv lJarvis & Jainl d2004l) . Yet, here we want to include 
the extra information encoded in the variation to constrain the 
galaxy field PSF models, and will therefore account for the non- 
linear dependence below. The mean stellar half-light radius in 
each exposure, is plotted as a function of the first principal 
component coefficient in Fig. IB. 71 showing that a fourth-order 
polynomial fit is capable to describe the full non-linear varia- 
tion. 

In order to obtain a well constrained model for all p a with 
high spatial resolution, we jointly fit all stars from all j < N ex- 
posures with a model 

1=0 c,=l 

separately for both chips, where P ac ^pi Cl (x,y) indicates a fifth- 
order polynomial in the corresponding rescaled x, y coordinates, 
and I = with co, max = 1 and yoj = 1 corresponds to the sub- 
tracted mean data vector, now modelled with high spatial res- 
olution. We aim to fit the few stars in the galaxy fields with as 
few parameters as reasonably possible. Due to the dominant role 
of focus changes we hence use only the first principal compo- 
nent in our analysis / max = 1, but include up to fourth-order terms 



(ci,max = 4) inyij. This takes out the non-linear distortion visible 
in Figs. IB. 6l and IB. 71 and hence the bulk of the variation in the 
second principal component. This combination yields a total of 
(1 +4)x21 = 105 coefficients per PSF quantity and chip, which 
are very well constrained from a total of 5 x 10 5 stars per chip. 

For illustration we plot the field-of-view dependence of the 
high-resolution DRZ ellipticity model measured for r g = 1 .4 pix- 
els in Fig. IB. 81 where the left panel shows the mean PSF elliptic- 
ity (Z = 0), while the right panel depicts the first singular vector 
{I — 1). Note the slight discontinuity of the mean PSF ellipticity 
between the chips, which is likely caused by small hei ght differ- 
ences between the C CDs as reported by Krist (2003). See also 
Rhod es et al. I d2007h who measure a stronger discontin uity in the 
TinyT im PSF model but not for stars in COSMOS, and lJee et all 
(2007) who notice it in the PSF size but not ellipticity variation. 

To obtain PSF models for our COSMOS stacks, we fit 

e v e 2> r h °f a ^ stars m tne sm gl e COSMOS COR exposures with 
the PCA model (IB. 1 11 1 to determine the first principal compo- 
nent coefficient y\j for this exposure. We then average the corre- 
sponding DRZ-image PSF models of all exposures contributing 
to a tile, taking thei r relative dither offsets an d rotations into ac- 
count, as detailed in Schrab back et al.l d2007l) . 

We plot the time dependence of the estimated coefficie nt y\ j 
for both the COSMOS and stellar field exposures in Fig. |B.9| 
Note that HST has been refocused at several occasions to com- 
pensate long-term shrinkage of the OTA, with one correction 
by +4.2 microns being applied during the time-span of the 
COSMOS observations on 2 004 D ec 2 j fj To ease the compari- 
son, all plots shown in Figs. lB.5l to lB.9l have been created using 
a single PCA model determined from all star fields. Yet, to ex- 
clude any possible influence of the refocusing, we actually use 
separate (but very similar) PCA models for the two epochs in 
our weak lensing analysis. 

While the fit (IB.1U captures ~ 97% of the total PSF varia- 
tion in the stellar fields and metric defined above, it is important 
to realize that further PSF variations beyond focus changes do 
actually occur. These are indicated by the higher principal com- 
ponents and the additional scatter beyond the curved distortion 

20 http : //www . stsci . edu/hst/observatory/f ocus/ 
mirrormoves. html 
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Fig. B.7. Mean stellar half-light radius r* as a function of the 
first principal component coefficient y\j for the 700 iai4 stellar 
field exposures. At yy ^ 1 the telescope is optimally focused. 
The curve shows the best fitting fourth-order polynomial fit. The 
outliers are caused by crowded fields with very broad stellar lo- 
cus. 

in the second principal component. The subdivision of fields ac- 
cording to the velocity aberration plate scale factor VAFACTOR in 
Fig. IB. 61 which depends on the angle between the pointing and 
the te lescope orbital velocity vector (see e.g. ICox & GiU iland 
l2002h . indicates that these distortions are not random but may 
be coherent for neighbouring fields observed under similar con- 
ditions. This is not surprising given that HST undergoes sub- 
stantial temperature changes and the relative angle towards the 
sun may lead to pointing-dependent effect^]. For a survey like 
COSMOS, where neighbouring fields have often been observed 
under similar conditions, we hence expect coherent residual PSF 
distortions beyond the one-parameter model introduced here. 

These residuals cannot be constrained reliably from the few 
stars present in a single ACS galaxy field, as one would have 
to fit > 10 principal components given the slow decline of the 
/ > 3 eigenvalues CFig. IB.5t . However, under the assumption that 
they are semi-stable for fields observed under similar conditions, 
we can constrain these PSF residuals by combining the stars of 
multiple COSMOS tiles taken closely in time. 

For this purpose, we split the COSMOS fields into 24 
epochs. Within each epoch the data were taken closely in time, 
with the same orientation and similar sun angles. The only 
exception are two tiles which were reobserved between 2005 
Oct 28 and 2005 Nov 24 due to previous guide-star failures 
dKoekemoer et al.l l2007). which we add to epochs observed one 
year earlier under similar conditions. Within each epoch we 
combine the stars of all tiles and compute residuals of the DRZ 
image PSF quantities e\, e* 2 , q\,q* 2 , T* by correcting their values 
measured from the stacks with the (dithered and averaged) mod- 
els ( IB. lll l. We then fit these residuals as 

p<:f(x,y) = P%(x,y) + —jP ( ^y) (B.12) 

where ~j has been averaged between the four exposures con- 
tributing to the stack /', and P <2 ' and P (2) , indicate second-order 
polynomials in the rescaled coordinates 5t, y determined for both 

21 Note that the actual impact of velocity aberration on object shapes 
is negligible for our analysis, as long as it is properly accounted for in 
the image stacking, as done by MultiDrizzle. 



chips together. Here we assume that the additional PSF varia- 
tions are in principle stable during each epoch, but their impact 
might depend on the actual focus position and hence y\]. Note 
that we do not use higher-order polynomials in x, y or non-linear 
powers of yfy, as we would otherwise risk over-fitting for epochs 
with few contributing exposures. Yet, we tested slightly higher 
orders for those epochs containing sufficiently many stars, yield- 
ing nearly unchanged results. In general we found that the fit- 
ted coherent PSF residuals are small, with a mean rms model 
ellipticity of 0.3% (for r g = 1.4 pixels). However, some epochs 
showed somewhat enhanced ellipticity residuals, with two exam- 
ples given in Fig. IB.10l motivating us to include this extra term in 
the galaxy PSF correction. In contrast we found that the residuals 
for T* are negligible. Also note that we deviate from our philoso- 
phy to obtain purely exposure-based models at this point, which 
is justified by the small and smooth (low-order) corrections ap- 
plied, which are only marginally affected by dithering. 

B.6. Galaxy correction and selection 

We measure galaxy shapes and correct for PSF effects as de- 
tailed in the previous subsections. We then select galaxies with 
cuts rh > l.2r^ mdx , where r^ mdx is the maximum half-light ra- 
dius of the 0.25 pixel wide, automatically determined stellar lo- 
cus in the image, S/N > 2.0, and Tr[P«]/2 > 0.1, identical to the 
cuts applied to the simulated data in App. IB.2l We also reject 
saturated stars and galaxies containing masked pixels (Sect.[2]i. 
In order to correct galaxy shapes for spurious CTI ellipticity, 
we fit the PSF anisotropy-corrected galaxy ellipticity component 
e\ m — e\ - P\™q*p with the power law model 

ca , ga] _ _ / FLUX \" F / r f r*/SKYr s 
Cl gl U0 3 e-/ ^ 3 pixel j ^0e-j 

with the mean sky level of the contributing exposures SKY, 
the time t = MJD - 52340 since the installation of ACS, the 
number of y-transfers ytrans, the SExtractor flux-radius rf 
(FLUX_RADIUS), and the mean integrated flux per exposure 
FLUX = 7^ • FLUX_AUT0 measured by SExtractor. We scale 
the latter with the mean exposure time per exposure as the 
stacks are in units of e~/ s . This model is similar to the one 
employed by Rhod es et al.l (120071) . but additionally accounts for 
the sky background-dependence of CTI effects and allows us 
to separate the dependence on galaxy flux and size. Despite 
the similarity to the stellar model dB.5b . we do not deter- 
mine a common CTI model for the typically bright stars and 
faint galaxies, as a simple power law fit is not guaranteed to 
work well over such a wide range in S/N. Considering all 
selected COSMOS galaxies we determine best fitting param- 
eters (e° v F,R,S) = (0.0230,0.134,0.638, 1.46). The correction 
for field distortion leads to a mean rotation of the original y-axis 
and hence readout direction in ACS stacks and DRZ exposures 

by <fi 2.5°. Thus, CTE degradation has also a minor effect on 

the <?2 ellipticity component, which we account for in both the 
galaxy and stellar correction as e" 6 = tan {2<p) ef e -0.087 ef e . 
Note that CTI affects an image after convolution with the PSF. 
Hence, one would ideally wish to correct for it first. Yet, in or- 
der to determine the impact of CTI, we need to correct for PSF 
anisotropy first, which would otherwise dominate the mean ei el- 
lipticity. We then subtract the CTI model (IB . 1 31 > and compute the 
fully corrected galaxy ellipticity ejf = (2/Tr[P«])(e™ - ef ,gal ) 
with ( IB.3l l. which is an unbiased estimate for the shear g a if 
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Fig. B.8. PCA PSF model ( IB.1U for the DRZ field-of-view ellipticity variation measured with r„ 
the mean ellipticity (I = 0), whereas the right panel depicts the first singular vector (Z = 1, q 
changes, for an arbitrary scale y\ = -3.0 (for positive y\ the ellipticities are rotated by 90°). 
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Fig. B.9. Temporal variation of the first principal component coefficient y\j, which is related to the HST focus position, measured in 
the stellar field and COSMOS exposures. The long-term shrinkage of the OTA is well visible as a decrease in the mean y\j, which 
was compensated with the marked focus adjustments. The large spread at a given date is not caused by measurement errors but 
orbital breathing leading to substantial short-term focus variations. 



©, © are taken into account. As it may be easier applica- 
ble for non-KSB methods, we also quote best-fitting parameters 
(e°,F,R,S) = (0.0342,0.068, 1.31, 1.26) if the actual shear esti- 
mates are fitted instead of the PSF anisotropy-corrected elliptici- 
ties, where the difference is caused by the PSF seeing correction 
blowing up the CTI ellipticity. 

As a test for residual instrumental signatures we create a 
stacked shear catalogue from all COSMOS tiles. Doing this, 
we marginally detect a very weak residual shear pattern, which 
changes with cuts on ~j. To quantify and model this residual pat- 
tern, we fit it from the PSF anisotropy and CTI-corrected galaxy 
ellipticities with a focus-dependent, second-order model (IB. 1 21 > 
jointly for all fields, yielding a very low rms ellipticity correction 



of ~ 0.003. One possible explanation for these residuals could 
be the limited capability of KSB+ to fully correct for a complex 
space-based PSF, despite the very good performance on the sim- 
ulated data in App. IB.2l Alternatively the limited number of stars 
per field may ultimately limit the possible PSF modelling accu- 
racy. In order to assess if these residuals have any significant 
impact on our results, we have performed our science analysis 
twice, once with and once without subtraction of this residual 
model. The resulting changes in our constraints on cr g are at the 
2% level, which is negligible compared to the statistical uncer- 
tainties. Also the E/B-mode decomposition is nearly unchanged 
(Fig.|4j». We only detect a significant influence for the star-galaxy 
cross-correlation, which is strictly consistent with zero only if 
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COSMOS DRZ: epoch 6 residual model r =1.4 COSMOS DRZ: epoch 16 residual model r =1.4 
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Fig. B. 10. Examples for the residual ellipticity model dB.12l > determined after subtraction of the 1 -parametric PCA model (1B.1U 
from the stellar ellipticities measured in COSMOS stacks with r„ — 1.4 pixels. The left {right) plot has been determined from all 
COSMOS fields with 732 < t < 735 (950 < t < 954.5), where t = MJD - 52340. Each whisker represents the residual ellipticity 
model for one star in the epoch. In some cases the model appears to be discontinuous due to the dependence on ~j or focus. Note 
the different scale compared to Fig. IB. 81 



this correction is applied, but even without correction it is negli- 
gible compared to the expected cosmological signal (Pig.|5j. 

As last step in the catalogue preparation, we create a joint 
mosaic shear catalogue from all fields, carefully rejecting dou- 
ble detections in neighbouring tiles, where we keep the detection 
with higher S IN and refine relative shifts between tiles. In the 
case of close galaxy pairs with separations < U.'5 we exclude the 
fainter component. Our filtered shear catalogue contains 472 991 
galaxies, corresponding to 80 galaxies/arcmin 2 , with a mean 
ellipticity dispersion per component cr e a = 0.31. To limit the 
redshift extrapolation in Sect. l2.2.2l we apply an additional cut 
z'gi4 < 26.7, leaving 446934 galaxies, or 76 galaxies/arcmin 2 . 

We rotate all shear estimates to common coordinates, and 
accordingly create a joint mosaic star catalogue for the analysis 
in Sect.@] 
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